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In  this  report  is  described  an  investigation  erf  the  i>tres2  distribu¬ 
tion  within  the  cross-sectica  of  e  gravity  dea.  The  analysis  was  performed 
by  c  agtccsetic  digital  cccputer,  using  the  finite  desert  szetbed .  The 
principal  purpose  of  the  study  was  to  detezaine  the  effect  on  the  stress 
distribution-  of  a  vertical  crack  extending  free  the  center  of  the  base 
through  cost  of  the  height  of  the  section-  The  dsn  was  sub  jeer  ad  ce 
dead  load,  live  load,  sad  tharsal  stress  conditions;  a  total  of  17 
different  loading  and  crack  ccafigiratioss  were  considered,  including 
cases  in  which  the  foundation  aaterial  was  assisted  to  be  oxthotrcpically 
clastic. 

Results  e£  the  studies  sre  presented  graphically  in  the  fois  of 
e tress  contours  and  boundary  displecexent  diagrams  for  each  case¬ 
in  addition,  figures  are  presented  depicting  the  s%>cst  of  opening  or 
dosing  of  the  crack  resulting  frees  the  various  load  conditions. 

Although  the  results  ahoy  that  definite  stress  concentrations  result 
frees  the  presence  of  the  crack,  the  uaxitars  stresses  that  develop  are 
not  severe,  being  less  than  about  5G0  psi  in  all  cases. 

Subsequent  to  the  completion  cf  the  original  report,  a  suppi enen t ary 
investigation  was  carried  out  in  which  additional  conditions  of  loading 
and  of  crack  configuration  were  considered.  Included  among  tnese  were 
cases  in  which  the  crack  was  assured  to  extend  through  the  upper  surface 
of  the  dac,  effectively  separating  it 


into  two  blocks.  Because  the 


ia 


filling  of  the  reservoir  behind  such  a  section  causes  the  crack  to  close 
gradually  from  the  top  downwards,  it  was  necessary  to  develop  an  in¬ 
cremental  analysis  procedure  for  the  full  crack  height  cases.  The  results 
of  this  supplementary  investigation  are  presented  here  as  a  supplement 
to  tne  origins^,  report  (see  p.  35).  They  shew  that,  although  the 
presence-  of  a  crack  extending  through  the  full  height  of  the  section 
leans  to  a  moderate  increase  of  stress  concentration  at  the  base,  the 
stresses  still  are  not  severe — only  about  600  psi  maximum. 
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TEE  STRESS  DISTRIBUTION  OF  NORFORK  DAM 


INTRODUCTION 


Origin  of  the  Problem 

The  distribution  of  stress  which  actually  exists  in  massive  concrete 

gravity  dams  has  long  been  the  subject  of  conjecture.  A  number  of 

(1,2,3)* 

experimental  studies  have  been  carried  out  *  *  both  in  the  laboratory 
and  in  the  field,  in  an  attempt  to  establish  the  order  of  magnitude  of 
the  errors  inherent  in  the  simplified  assumptions  which  form  the  basis 
for  design.  The  results  of  these  investigations  have  shed  considerable 
light  on  the  problem  and  bear  out  the  adequacy  of  the  design  assumptions 
in  normal  cases.  However,  little  is  known  about  the  distribution  of 
stress  in  dsuns  in  which  conditions  differ  radically  from  the  normal. 

One  condition  which  is  known  to  exist  in  a  number  of  structures, 
but  the  effect  of  which  has  not  been  studied  previously,  is  the  presence 
of  a  vertical  crack  or  system  of  cracks,  extending  in  planes  parallel 
to  the  longitudinal  axis  of  the  dam.  Such  a  crack,  if  it  extended 
through  a  major  portion  of  the  height  of  the  dam,  obviously  could  have 
a  significant  effect  on  its  structural  action  in  resisting  water  loads, 
and  conceivably  could  reduce  materially  its  overall  factor  of  safety. 

The  purpose  of  the  investigation  described  in  this  report  was  to 
determine  the  influence  of  such  a  vertical  crack  on  the  stress  distribu¬ 
tion  and  deformations  of  an  actual  concrete  gravity  dam.  The  subject  of 
the  study  was  the  Norfork  Dam,  a  230  foot  high  structure  built  on  the 
North  Fork  River  in  Arkansas  by  the  United  States  Army  Corps  of  Engineers. 
Extensive  cracking  was  noticed  in  this  structure  during  its  construction. 


*Numbers  in  parentheses  refer  to  similarly  numbered  References  listed  at 
the  end  of  this  report. 
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and  a  rather  extensive  survey  of  the  distribution  and  else  of  the  principal 
cracks  has  been  carried  out  during  tha  18  yaars  since  its  completion,  (hi 
the  basis  of  this  study  it  was  apparent  that  Monolith  1$  vat  ona  of  tha 
siost  savarely  cracked  portions  of  tha  structure  and  it  was  selected  as  c 
subject  for  mathsaatical  analysis. 

Within  this  monolith  was  one  principal  crack,  located  approximately 
in  mid-width  and  extending  vertically  ever  most  of  its  height.  The 
specific  objective  of  this  study  was  to  determine  haw  such  a  single  crack 
in  an  otherwise  solid  section  would  affect  the  distribution  of  stress 

i 

within  the  section.  The  vertical  extent  of  the  crack  could  not  be 
established  with  great  reliability,  thus  this  factor  was  selected  as  one 
of  the  prindpsl  variables  of  the  study. 

Analytical  Procedure 

A  second,  and  also  very  valuable,  purpose  of  the  investigation  was 

(A) 

the  evaluation  of  a  recently  developed  method  of  plane  stress  anelysis 
This  method,  which  is  known  as  the  "finite  element  method",  takes  ad¬ 
vantage  of  the  tremendous  computational  capacity  of  modern  high-speed 
digital  computers.  It  had  been  tasted  previously  in  relatively  minor 
acadwlc  problems,  but  its  effectiveness  in  solving  c  full  scale  practical 
problem  had  not  yet  been  demonstrated. 

The  fundamental  concept  of  the  finite  element  method  was  developed 

(5) 

originally  in  an  attempt  to  extend  and  generalize  the  "lattice  analogy" 

Which  had  been  proposed  as  a  method  of  plane  stress  analysis  by 

Hrennikoff^aad  McHenry^  over  twenty  years  ago.  An  interesting  aide 

issue  of  this  investigation  was  the  offer  of  Hr.  McHenry *  now  working 

with  the  Portland  Cement  Association,  to  make  a  similar  analysis  by  the 

Director  of  Development,  Portland  Cement  Association  Research  and  Development 
Laboratories ,  Skokie,  Illinois. 


lattice  analogy  procedure  (which  had  not  previously  been  programmed  for 
the  digital  computer).  This  analysis  provided  an  independent  check  on 
the  results  obtained  by  the  finite  element  method,  and  also  made  possible 
a  comparison  of  the  relative  merits  of  the  two  procedures. 

Scope  of  the  Investigation 

The  geometry  of  the  cross°section  considered  in  the  investigation 
and  the  properties  of  the  concrete  and  foundation  rock  were  specified  by 
the  Contracting  Office.  Loadings  which  were  considered  included  dead 
load,  live  load  (hydrostatic  pressure),  and  the  effects  of  temperature 
changes.  Stress  distributions  and  deformations  were  calculated  due  to 
these  load  conditions  applied  in  various  combinations,  as  agreed  upon  by 
the  investigators  and  the  Contracting  Office.  The  initial  width  of  crack 
opening  was  defined  under  a  combined  temperature  change  and  dead  load 
condition.  The  additional  influence  of  the  live  load  in  closing  the 
crack  and  stressing  the  structure  was  then  evaluated.  Because  the 
foundation  rock  appeared  to  have  a  considerable  amount  of  Jointing  in 
horizontal  planes,  it  seemed  probable  that  its  effective  modulus  of 
elasticity  in  the  vertical  direction  was  considerably  less  than  in  the 
horizontal  direction.  Therefore,  the  effects  of  an  orthotropic  foundation 
was  studied  in  two  cases. 

All  of  the  analyses  were  performed  on  the  IBM  704  electronic  digital 
computer  operated  by  the  University  of  California  Computer  Center.  The 
principal  computer  program  used  in  the  analyses  was  developed  by  the 
investigators,  prior  to  the  initiation  of  this  contract,  on  a  research 
grant  from  the  National  Scieilce  Foundation.  However,  a  few  minor  program 
modifications  and  extensions  were  made  during  the  course  of  this  study. 

Stress  distribution  and  displacement  results  were  printed  by  the  computer  in  tab¬ 
ular  form,  but  were  subsequently  portrayed  graphically  for  easier  comprehension. 
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METHOD  0?  AHALYSIS 

The  basic  problem  to  be  solved  in  this  investigation  was  the 
determination  of  the  stresses  and  displacements  in  a  plane-stress  system 
of  arbitrary  dimensions  and  properties;  due  to  the  application  of  tempera¬ 
ture  changes,  external  loads,  and  its  own  dead  weight.  The  finite  element 

(4) 

method  which  was  employed  in  the  analysis  has  been  described  previously 
so  the  basic  principles  of  the  procedure  will  be  reviewed  here  only  briefly. 
However.,  several  aspects  of  the  present  problem  required  the  development 
of  new  concepts  and  techniques  which  will  be  described  in  full. 

The  Finite  Element  Approximation 

The  basic  approximation  introduced  in  this  method  of  analysis 
is  the  representation  of  the  actual  plane-stress  system  as  an  assemblage 
of  triangular  plate  elements,  interconnected  only  at  the  corners  or 
nodal  points,  and  each  subjected  to  a  state  of  plane  stress.  This 
idealization  replaces  the  actual  elasticity  problem  by  a  structural 
analysis  problem.  The  advantage  of  the  substitution  is  that  the  conditions 
of  compatibility  and  equilibrium,  which  must  be  represented  by  partial 
differential  equations  in  the  elasticity  problem,  may  be  expressed  by  a 
set  of  simultaneous  algebraic  equations  in  the  structural  problem.  This 
problem  may  then  be  solved  by  standard  methods  of  structural  analysis. 

A  very  coarse  mesh  (and  therefore  quite  crude)  idealization  of  a 
dam  cross-section  is  shown  in  Fig.  1.  The  triangular  elements  are  joined 
only  at  their  corners  and  the  loads  are  also  applied  only  at  these  nodal 
points.  Deformations  of  the  assemblage  depend,  of  course,  not  only  on 
the  stiffnesses  of  the  individual  elements,  but  also  on  the  manner  in 
which  they  are  coupled.  Thus  it  would  appear  that  the  idealization 
would  be  excessively  flexible,  due  to  the  fact  that  continuity  between 


the  elements  is  imposed  at  so  few  points.  However.,  this  incomplete 
continuity  condition  is  compensated  by  the  assumption  that  the  stresses 
within  each  element  must  be  uniform:  i.e.,  the  normal  stresses  in  the 
horizontal  (  o*  )  and  vert,  .a!  '&y)  directions,  and  shearing  stress  ( '£  ) 
are  constant  within  each  element.  As  a  result  of  this  assumption,  the; 
edges  of  each  element  remain  straight  under  load,  and  therefore  continuity 
is  achieved  not  only  at  the  nodal  points  but  also  along  the  element 
boundaries  as  well. 

It  should  be  noted  that  one  of  the  major  advantages  of  the  finite 
element  approximation  is  that  the  properties  of  each  element  may  be 
selected  independently.  Thus  different  moduli  of  elasticity,  Poisson's 
ratios,  unit  weights,  thermal  changes,  etc.,  can  be  imposed  in  each 
element  without  increasing  the  complexity  of  the  analysis.  This 
characteristic  of  the  method  is  particularly  advantageous  in  the  analysis 
of  dams  where  a  definite  discontinuity  of  properties  is  likely  to  exist 
between  the  foundation  rock  and  the  concrete  of  the  structure. 

The  Matrix  Analysis  Procedure 

The  finite  element  idealization  of  a  plane  stress  problem  constitutes 
a  highly  interconnected  structural  system.  For  such  systems,  the  displace¬ 
ment  method  of  structural  analysis  has  proven  to  be  most  effective;  and 
for  large  systems,  the  use  of  matrix  algebra  leads  to  the  most  efficient 
computational  processes.  A  general  description  of  the  displacement  method 

formulated  as  a  sequence  of  matrix  operations  has  been  published  previous- 

,  (8) 

ly  ,  but  the  principles  will  be  reviewed  here  for  the  convenience  of 
the  reader. 

The  analysis  begins  with  a  description  of  the  structure  and  its 
loading  expressed  in  matrix  form.  The  properties  of  the  structure  are 
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defined  in  two  matrices;  one,  called  the  element  stiffness  matrix,  expresses 
the  elastic  properties  of  each  of  the  triangular  elements  making  up  the 
structure;  the  second,  called  the  displacement  transformation  matrix,  de¬ 
fines  the  geometry  of  the  assemblage.  In  addition,  the  listing  of  the 
vertical  and  horisontal  load  components  applied  to  each  nodal  point  comprises 
the  load  matrix.  These  three  matrices,  which  will  be  designated  fkj  , 
f  a j  ,  and  |r|  respectively,  are  discussed  further  m  the  following  sections 
of  this  chapter. 

When  the  properties  of  the  structure  have  been  expressed  in  matrix 
fora,  the  analysis  proceeds  as  follows.  First,  the  stiffness  of  the  assembled 
Structural  system  [kJ  is  obtained  by  multiplying  together  the  element  stiff¬ 
ness  and  displacement  transformation  m  trices: 

CkI-WTMM  « 

in  which  f_a]T  is  the  transpose  of  [aj  .  The  structural  stiffness  matrix 
[k]  defines  the  nodal  force  components  in  the  assembled  structure  which 
result  from  unit  displacements  of  the  nodal  points.  Representing  these 
elastic  force  components  and  nodal  displacements  by  the  column  matrices 
W  and  [r]  ,  respectively,  this  force-displacement  relationship  may 

be  expressed: 

tel  =  MM  c» 

However,  equilibrium  considerations  require  that  the  elastic  forces  {q] 
developed  at  the  nodal  points  must  equal  the  externally  applied  loade 

.  Thus  Eq.  2  may  be  rewritten  in  terms  of  the  external  load  matrix: 

..  tel  =  MM 


(3) 
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In  a  typical  problem  the  loads  {r]  are  given,  while  the  displace- 
meats  £rj  are  to  be  determined,  thus  the  next  step  in  the  analysis  is 
the  inversion  of  the  stiffness  matrix  |k]  to  obtain  the  flexibility  of 
the  structure  [f]  : 

Lf]  =  [k]"  (4) 

Multiplying  both  aides  of  Eq.  3  by  the  flexibility  matrix  then  yields 
an  expression  for  the  displacements  in  terms  of  the  loads, 

M  »  !/]{*}  <5> 

Eq.  5  thus  provides  one  part  of  the  final  analytical  results  which  were 
required:  the  deformations  of  the  structure. 

After  the  nodal  displacements  have  been  determined t  the  deformations 
of  the  individual  plate  elements:  denoted  by  the  column  matrix  ^U"} 
can  be  determined  by  applying  the  displacement  transformation  matrix, 

M  *  MW  <6> 

Multiplying  these  element  deformations  by  the  element  stiffness  matrix 
then  yields  the  nodal  point  forces  {s}  in  each  element, 

.  W  c  0]W  c  WWW  (7> 

(It  may  be  noted  that  Eqs.  6  and  7  essentially  constitute  definitions 
of  the  matrixes  [a]  and  [k]  ). 

In  the. analysis  of  a  structural  system  which  is  composed  of  bar  type 
elements,  the  nodal  point  forces  would  represent  the  final  results.  How¬ 
ever,  in  the  plane  stress  problem  considered  here,  the  distributed  stresses 
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within  the  elements ,  {o'}  ,  rather  than  the  resultant  nodal  forces,  are 
of  primary  interest.  Thus  one  final  step  is  needed,  the  transformation 
from  nodal  forces  to  element  stresses.  Designating  the  stress  transfor¬ 
mation  matrix  by  the  symbol  [w]  t  this  operation  is  performed  as  follows 

to-}  =  [m]{s}  =  H(kJ  {tr}  (8) 

Finally,  representing  the  above  matrix  product  by  a  single  symbol 

[m]  2  HW  <9> 


Eq.  8  may  be  written 


H  •  MWid  <w) 

in  which  [M]  performs  the  function  of  transforming  from  element 
deformations  to  element  stresses.  The  entire  analytic  procedure  is  now 
contained  in  Eqg.  1,  4,  5,  and  10, 

Although  the  procedure  described  above  is  straight-forward,  and 
provides  a  convenient  method  for  the  analysis  of  small  to  moderate-size 
systems,  the  inversion  of  the  stiffness  matrix  (Eq,  4)  is  a  tremendous 
task  when  large  systems  (involving  over  perhaps  100  nodal  points)  are  to 
be  analyzed.  In  order  to  avoid  this  problem,  an  iterative  solution  of 
Eq,  3  for  the  nodal  point  displacements  was  substituted  for  the  inversion 
process.  Using  this  approach  the  procedure  is  capable  of  treating  systems 
having  as  many  as  600  nodal  points. 
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Element  Stiffnesses 

Isotropic  Elasticity  -  The  method  used  in  deriving  expressions  for  the 
stiffnesses  of  triangular  plate  elements  composed  of  Isotropic  elastic 

(4) 

material  has  been  described  previously  so  only  the  results  will  be 
presented  here.  The  geometry  of  a  typical  element  is  defined  in  Fig.  2. 


FIG.  2  TRIANGULAR  PLATE  ELEMENT  DIMENSIONS 
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Its  scifn«s  satrin  is  presssted  is  Sq,.  11,  b^lcs: 


Sc _ 

20-^X^-feo.j 


fcV'fdk 


(Sfasta'kd)  ! 


(11) 


ic  vhich 

E  =  Modules  of  Elssciciey 
t)  ;  ratios 


In  the  preceding  diecsssien  c£  the  fisita  element  cached,  ic  has 
been  assumed  chst  the  system  co  be  solved  is  subjected  cc  a  co edition 
of  piece  stress,  ard  2q-  11  applies  co  this  condition.  to  the  other 
hand,  there  say  be  scce  jus tificacico  for  as  seal  eg  chat  e  srcritT  das t  is 
in  a  state  of  plane  a  train-  Ifce  plane  stress  stiffness  matrix  of  3q.  11 
can  be  applied  in  a  plane  strain  analysis  if  modifies?  elastic  properties 
are  introduced  in  the  equations  as  follows  • 


fc— 

1)* 


(12) 


In  these  equations,  the  "starred'’  valuer  are  the  modified  properties 
to  be  used  in  applying  Eq.  11  to  a  plane  strain  analysis,  vhiie  the 
uns tarred  values  are  the  actual  properties. 

In  the  present  investigation,  the  plane  stress  equations  vere  used 


X  J 

excicst-aaly  btianse  1 t  was  ssrsxed  chat  the  astfe od  of  c%»tpn:lo,5 
lie  large  oocolitrs)  *<7-ld  lead  sere  closely  co  suet  c radi cions .  Hxe-et 
it  nay  also  be  coted  that  ter  the  values  of  Poisson1  s  ratio  :scO  -r  tils 
2 tody  {  i)  =  G .  1 7 }  the  difference  between  the  -o.o3l'l ee  and  art ^e'.  p r ope: 
ties  (  E*  ~  1.03  E  :  ~0*  ° i  20l)  *  Is  well  within  the  ..rite  of 
accuracy  vita  which  these  values  caa  be  estimated  Thus  'he  quean-ai  as 
Co  whether  a  condition  of  plane  stress  or  plane  strain  eaists  la  only 
of  academe  interest. 

Orthotrooic  elasticity  -  In  order  to  establish  the  Influence  of  an 
orthccropic  fo-usdsticra  aaterlai  on  the  stress  distribution  in  the  dan,  it 
was  necessary  to  develop  the  stiffness  satrix  for  ea  orthotropic  triangular 
eleaent.  The.  derivation  of  the  orthotropic  stiffness  sea  cm  followed 
exactly  the  sane  procedure  used  ta  deriving  the  isotropic  stiffness  matrix 
except  that  the  stress- strain  relationship  for  the  orthotropic  material 
was  substituted  for  the  isotropic  stress-strain  relationship  used  in 
Heference  &.  Thus  the  first  requirement  was  the  derivation  of  the  ortho- 
tropic  stress-strain  relationship. 

For  this  purpose,  it  was  assumed  that  the  orthotropic  material  wag 
actually  a  layered  system  consisting  of  two  isotropic  materials  of 
different  stiffnesses.  and  E2.  as  shown  in  Figure  3. 


IIG.  i  ASSUMED  ORIfc 3TR0FIC  MATERIAL 

The  softer  material  (E„)  was  assigned  to  occupy  only  a  small  proportion 
r,  of  the  total  thickness,  and  It  was  further  assumed  that  the  ratios 
o'  the  elastic  moduli  of  the  two  materials  was  the  same  as  the  ratios 
of  their  Poisso n' a  ratios 


t>r 


(13) 


The  effective  elastic  properties  of  this  layered  material  were,  determined 

by  applying  successively  stresses  CT  C T  ■■  and  7T  .  and  computing 

x‘  y  xy' 

the  resulting  strains  6  6  ,  and  S'  . 

x  y  xy 

Applying  the  horizontal  stress^  (T*  the  effective  horizontal 

X  ’ 


modulus  is  found  to  be 


<rK 


where 


U4) 


while  the  effective  Poisson's  ratio  is  given  by 


A 


-  i ),n-r*) 


(15) 


SLaila.-Ly,  applying  the  vertical  stress.  C£  ,  the  effective  vertical 
modulus  is 


% 


_  Ez 

£u.  £2  .  ✓-  * 

*  -p  ^  • 

w  * 


(16) 


and  the  corresponding  Poisson  3  rario  is 


1)  =  -  _<L 

V*A  ~  F. 

ch 


=  l)*-! 

St  > 


(17) 


where 


/ 


Finally,  applying  the  shear  stres3;,  s  the  effective  shear  modulus 

is  found  to  be 

= 

2(ju.  + 


<5 


(18) 


These  elastic  properties  may  be  summarized  by  the  ortho-tropic. 


stress-strain  relationship,  which  follows- 
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In  the  derivation  of  the  element  stiffness.,  the  inverse  iorm  of  this 
stresa-s train  relationship  :.a  required  as  follows 
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Employing  this  orthotropic  stress- strait-  re.1  at  ionshtp  in  the  procer.ce 
described  in  Reference  a  for  the  derivation  of  the  tri-a'igulii  eiem-st-' 
stiffness  leads  to  the  folio- log  final  result  for  the  orthotopic  plate 
element  stiffness  matrix.  L&-] 


£t 


/>b 1^4 
y^Ay^A 

%«A+Za*bj 


a if  *^bjL 


(21 ) 


_  p  -  A 
?  ~  2.(//x+'0*) 


where 
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Complete  Element  Stiffness  Matrix  -  The  complete  elment  stiffness  matrix 
is  formed  merely  by  arranging  the  individual  element  stiffnesses  as 
submatrices  on  the  diagonal  of  the  complete  matrix.  Thus,  if  the  individual 
element  stiffness  matrices  of  elements  a,  b,  c,  .  .  .  n  are  designated 
ka,  kc»  •  •  •  kn  £ke  complete  matrix  may  be  written 


w 


(22) 


Element  Stress  Matrix  -  It  was  noted  above  that  the  stresses  in  an 
element  are  of  greater  interest  than  the  nodal  forces,  and  that  the 
matrix  [m]  would  serve  to  transform  from  the  element  deformations  to 
the  stresses.  This  stress  transformation  matrix,  which  was  not  presented 
in  the  previous  publication,  is  given  below  for  any  arbitrary  element  "i". 


[Mil  a 


_ E _ 


bi;  -  b,  t>a; 
-1%  0.; 
-§  bj 


where 
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(23) 


may 


The  transformation  of  stress  for  all  elements  in  the  structure  thus 
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ba  accomplished  by  means  of  a  diagonal  matrix  made  up  of  the  individual 
element  stress  transformations, 


Mb 


M* 


(24) 


Although  the  normal  and  shear  stresses  (  <3^  ,  Otj,  £  )  completely 

define  the  state  of  stress  in  the  elements,  it  frequently  is  of  interest 
to  knew  the  principal  stresses  and  Cfj.  and  their  directions  @  . 

In  this  investigation  these  principal  stress  values  were  determined 
from  the  stress  values  related  to  the  x,  y  coordinate  axes  by  the 
standard  transformation  formulas  of  elementary  mechanics. 

Displacement  Transformation  Matrix 

The  displacement  transformation  matrix  serves  to  relate  the 
deformations  of  the  finite  elements  to  the  displacements  of  the  nodal 
points,  as  is  shown  by  Eq.  6.  The  derivation  of  this  matrix  will  be 
discussed  with  reference  to  Fig.  4,  which  shows  a  portion  of  a  triangular 
element  system.  The  form  of  the  displacement  transformation  can  be 
demonstrated  by  establishing  the  rationship  between  the  displacements 
of  nodal  points  I,  J,  and  K,  and  the  resulting  deformations  of  element 
**n”.  Similar  relationships  could  be  established  between  any  nodal 
displacements  and  the  deformations  of  the  elements  connected  thereto. 
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b.  NODAL  DISPLACEMENTS  DUE  TO  RIGID  BODY  ROTATION 


FIG.  4 


TYPICAL  DISPLACEMENT  TRANSFORMATION  ANALYSIS 
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As  an  example,  a  unit  vertical  displacement  at  joint  J  is  indicated 
in  Fig.  4a,  anc  tha  resulting  distorted  configuration  of  element  MnM  is 
shown  by  the  dashed  lines.  This  distortion  pattern  violates  the  support 
conditions  which  were  assumed  for  the  element  in  deriving  its  stiffness 
(see  Fig.  2).  To  return  the  elenent  to  the  required  support  conditions, 

*  gigld  body  rotation  is  applied,  as  shown  in  Fig.  4b.  The  actual 
distortion  shown  in  Fig.  4a,  can  now  be  expressed  in  terms  of  the  required 
support  conditions  as  is  indicated  in  the  second  column  of  the  matrix  below 


(25a) 


The  first  column  of  this  matrix  represents  the  relationship  between 
the  horizontal  displacements  of  joint  J  and  the  element  deformations. 
Since  this  displacement  does  not  violate  the  required  support  conditions 
for  the  element,  it  is  clear  that  U-i  *  directly. 

Eq.  25a  may  be  written  symbolically  as  follows 

M 


Similar  relationships  can  be  developed  for  the  displacements  of  nodal 
points  X  and  K  as  follows 


or  in  symbolic  £oss 


=  \  {^1  *■  L<,nK]{r*c} 


Eqs.  25a  and  25b  may  also  be  combined  Into  a  single  relationship t 

I 

expressed  symbolically  as  follows 


1  <J"n  \  c  (flwr 


(26) 


Relationships  similar  to  Eq.  26  can  be  established  between  all 
uodal  point  displacements  and  element  deformations.  For  example,,  for 
the  simple  system  shown  in  Fig.  5,  the  complete  displacement  trans¬ 
formation  would  be  of  the  form 


'0* 

m  m 

Gilt,  Q&Z  Q 
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O  dj*;  o 
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O  «3tt  Qtn 
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O  0  «4K 

(27) 


in  which  each  of  the  symbols  is  as  defined  in  Eqs.  25a  and  25b „  with 
dimensions  selected  as  appropriate  to  the  particular  element. 
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Load  Matrix 

The  load  tutrix  la  aerely  a  listing  {r\  of  the  load  coaponent* 
applied  to  the  nodal  points  of  the  plats  alsaent  assemblage,  presented 
in  the  8«e  sequence  as  were  the  displacement  cccponents  £rj  in  the 
develop® ent  of  the  displaces ent  transformation  matrix.  Three  typer,  of 
loadings  were  considered  in  this  investigation,  dead  loads,  live  loads, 
and  thermal  loads,  each  of  which  will  be  discussed  separately  below. 

Bead  Loads  *  The  dead  weight  of  each  plate  element  is  given,  of  course, 
by  the  product  of  its  axe a,  its  thickness  (which  was  taken  as  one  foot) 
and  its  unit  weight,  ^  .  One  third  of  the  total  eleaent  weight  was 

applied  at  each  of  the  three  nodal  points  to  which  it  was  attached.  That 
the  total  dead  load  at  any  nodal  point  was  taken  ss  one-third  of  the 
weights  of  all  elements  attached  at  that  point,  applied  in  a  downward,  or 
negative  nyn  direction. 

Live  Loads  -  7m  calculating  the  hydrostatic  water  load  acting  at  each 
nodal  point,  it  was  assumed  that  the  pressure  was  distributed  against 
an  auxiliary  stringer  system,  as  shown  in  Fig.  6b.  The  reactions  at  the 
ends  of  the  stringers,  resolved  into  vertical  and  horizontal  conponsnts 
as  shewn  in  Fig.  6c,  then  constituted  the  live  loads.  For  any  given 
nodal  point,  the  live  load  depends  not  only  on  the  depth  of  water  at 
that  level,  but  also  upon  the  depths  at  the  points  above  and  below.  Where 
the  slopes  of  successive  element  faces  are  the  same,  as  is  the  case  in 
Fig.  6,  the  live  load  components  at  an  arbitrary  nodal  point  N  are  given  by 


F!G.  6  REPRESENTATION  OF  HYDROSTATIC  PRESSURES  BY  NODAL  FORCES 
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__  ^ 

~  ~~q~  ( dc-  d,, -*■  dU)6i»-<£-») 

I^JK  R'-t  -  Ccf  4* 

la  which  T*  ia  tJa  rail  weight  ef  water  a ad  the  other  symbols  sre 
defined  in  Fig.  6a. 

Thermal  leads  -  The  thermal  leads  stay  be  calculated  by  first  cal  cal  at  lag 
tbe  stresses  which  wedd  exist  in  the  detests  If  strains  due  to  tempera¬ 
ture  charges  were  fully  constrained.  In  a  plane  stress  system*  these 
stresses  are  glirea  by 


^  -  0^  =  <r  -  (29) 

in  which  cd.  r  thermal  coefficient  of  experts  ion 

AT  =  change  of  temperature 

The  nodal  forces  required  to  maintain  these  internal  stresses  in  each 
eletfflt  nay  then  oe  found  oy  simple  statics.  While  these  constraining 
forces  are  applied,  no  displacements  or  deformations  of  the  structure 
nay  occur,  and  the  stresses  are  given  directly  by  Eq.  29.  However,  the 
constraints  are  cot  really  present,  so  their  effect  trust  be  eliminated 
**y  applying  ecuai  and  opposite  nodal  forces.  These  reversed  constraining 
forces  are  the  thermal  Icaas  for  which  the  structure  is  analyzed. 


*In  a  plane  strain  system,  the  corresponding  expression  is 

__  _  E  cc  # 
c  O+^Xi-20) AT 

but  for  the  relatively  lew  Poisson5 s  ratio  assured  in  this  investigation, 
the  difference  between  this  and  Eq.  29  is  less  than  5  per  cent.  As  was 
mentioned  previously,  the  plane  stress  condition  was  assumed  throughout 
this  investigation. 
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For  the  element  shown  in  Fig.  2, 


they  would  be  expressed  by 


«c 

Rat 

i 

-K 

*y 

~  a 

(30) 


in  which  &  is  given  by  Eq.  29. 

The  displacements  resulting  from  these  thermal  loads  are  the 
true  thermal  displacements  in  the  system.  However,  the  stresses  resulting 
from  these  loads  are  only  the  part  of  the  thermal  stress  associated  with 
the  nodal  displacements.  The  total  thermal  stress  is  the  sum  of  the 
deformation  stress  resulting  from  the  thermal  loads,  plus  the  constrained 
stress  given  by  Sq.  29. 

Boundary  Displacements 

In  the  preceding  discussion  of  the  general  matrix  analysis  procedure, 
it  was  indicated  that  the  typical  analysis  requires  the  determination  of 
nodal  displacements  resulting  from  a  given  3et  of  applied  loads.  In  some 
cases,  however,  a  "mixed"  problem  may  exist  in  which  forces  are  specified 
at  only  part  of  the  nodal  points  and  certain  specific  displacements  are 
imposed  at  other  points.  Such  mixed  problems  may  be  solved  with  no  greater 
difficulty  than  cases  in  which  nodal  point  loads  are  the  only  stress  pro¬ 
ducing  effects. 

The  analysis  is  formulated  by  partitioning  the  nodal  point  forces 
into  two  submatrices:  the  specified  loads  and  the  unknown 

force  components  &}  (which  correspond  with  specified  nodal  displace* 
ments),  At  the  same  time,  the  nodal  point  displacements  {  r}  are 
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partitioned  into  submatrices  representing  the  specified  displacements 
^r^  and  the  unknown  displacements  (which  correspond  with  the 

specified  nodal  forces).  In  this  partitioned  form,  Eq.  3  may  be  written 


(31) 


in  which  the  stiffness  matrix  £kJ  has  been  partitioned  to  correspond 
with  the  partitioning  of  the  load  and  displacement  matrices.  Then  a 
formal  solution  for  the  unknown  displacements  is  obtained  from  the 

first  of  these  6ubmatrix  equations: 


W  =■  +  <«*> 

l.e.  {V\  =•  [K0„]"{RS$  -  [lC„y'[K.,](rs5  (32) 

Finally  the  unknown  forces  {iTj  may  be  determined  from  the  remaining 
submatrix  equation 

{X\  *  [K,0j[U$  +  (3lb) 

or  introducing  the  value  of  {0”|  from  Eq.  32, 

lX\  -  [Ki0][Keo]  {RS^[LK„]-[l<ij[KcJ  (33) 


Because  the  matrix  inversions  required  in  the  formal  solution  of 
a  large  structural  system  (the  calculation  of  Q^oo]  '  *-n  Eq.  33) 


is  a  formidable  problem,  an  iterative  solution  for  the  displacements  was 
employed  in  this  investigation,  as  has  been  mentioned  earlier.  The 
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solution  of  a  mixed  problem  in  which  certain  displacements  |>rs5  are  specified, 
does  not  require  any  modification  of  the  iterative  process.  No  iteration  is 
required  at  these  points  because  a  reaction  is  available  to  support  any 
unbalanced  forces  which  develop. 


29 


COMPUTER  PROGRAM 

The  digital  computer  program  performs  three  major  tasks  in  the 
complete  analysis  of  a  plane  stress  system.  First  the  stiffness,  displace- 
rnent  transformation,  and  load  matrices  are  formed  from  a  basic  numerical 
description  of  the  system.  Second,  Equation  3  is  solved  for  the  displace¬ 
ments  of  the  nodal  points  by  an  iteration  procedure.  Third,  the  internal 
element  stresses  are  determined  from  these  displacements.  The  general 
features  of  the  program  which  carries  out  these  operations  will  be  described 
in  this  chapter. 

Numerical  Procedure 

Before  presenting  the  sequence  of  operations  that  is  performed  by 
the  computer  program,  it  is  necessary  to  discuss  in  soma  detail  the  actual 
numerical  method  that  was  employed.  This  method  is  a  modification  of  the 
well-known  Gauss-Seidel  iteration  procedure  which,  when  applied  to 
Equation  3,  involves  the  repeated  calculation  of  new  displacements  from 
the  equation 


C«  1,0-1 


knc  n 


<»] 


(34) 


where  n  is  the  number  of  the  unknown  and  s  i~  the  cycle  of  iteration. 

The  only  modification  of  the  procedure  introduced  in  this  analysis 
is  the  application  of  Eq.  34  simultaneously  to  both  components  of 
displacement  at  each  nodal  point.  Therefore,  rn  and  Rq  become  vectors 
with  x  and  y  components,  and  the  stiffness  coefficients  may  be  expressed 
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in  the  form 


on 


(35) 


in  which  1  and  m  are  nodal  point  numbers. 

Over-Relaxation  Factor  -  The  rate  of  convergence  of  the  Ganss-Seidel 

procedure  can  be  greatly  increased  by  the  use  of  an  over-relaxation 
(9) 

factor  .  However,  in  order  to  apply  this  factor  it  is  first  necessary 
to  calculate  the  change  in  the  displacement  of  nodal  point  n  between 
cycles  of  iteration: 


(36) 


The  substitution  of  Equation  34  into  Equation  36  yields  for  the  change 
in  displacement 
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(37) 


The  new  displacement  of  nodal  point  n  is  then  determined  from  the 
following  equation: 


(vh) 


(38) 


where  is  the  over-relaxation  factor.  For  the  structure  considered  in 
this  report  it  was  found  that  a  value  of  ^  equal  to  1.86  gave  the  most 
rapid  convergence. 
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Physical  Interpretation  of  Method  -  There  is  important  physical 
significance  in  the  terms  of  Eq.  37.  The  term  (km)-l  is  the 
flexibility  of  nodal  point  n.  This  represents  the  nodal  point 
displacements  resulting  from  unit  nodal  point  forces,  and  can  be 
written  in  the  form  of  a  sub-matrix. 
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The  summation  terms  in  Eq.  37  represent  the  elastic  forces  acting 
at  nodal  point  n  due  to  the  deformations  of  the  plate  elements: 
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The  difference  between  these  elastic  forces  and  the  applied 
loads  is  the  total  unbalanced  force  which  in  sub-matrix  form  may 
be  written: 
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Eq.  38,  which  gives  the  new  displacement  of  nodal  point  n,  may 
now  be  rewritten  in  the  following  sub-matrix  form: 
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It  is  important  to  note  that  any  desired  nodal  point  displacement  r^ 
may  be  assumed  for  the  first  cycla  of  iteration.  A  good  choice  of  these 
displacements  will  greatly  speed  the  convergence  of  the  solution.  In 
fact,  if  all  displacements  were  assumed  correctly,  the  unbalanced 
forces  given  by  Sq.  41  would  be  zero  and  no  iteration  would  be  necessary. 
However,  in  a  practical  case  there  always  will  be  unbalanced  forces  in 
the  system  at  first,  and  the  iteration  process  continually  reduces  them 
toward  zero. 

Boundary  Condition  -  Eq.  42  is  valid  for  all  nodal  points  within  the 
structure;  however,  in  order  for  it  to  be  applied  to  boundary  nodal  points 
the  flexibility  coefficients  must  be  modified  to  account  for  the  specific 
types  of  restraint  which  may  exist.  Since  these  flexibility  coefficients 
are  independent  of  the  cycle  of  iteration,  this  modification  is  performed 
before  the  start  of  iteration.  There  are  three  possible  types  of  boun¬ 
dary  nodal  points  whose  flexibility  coefficients  must  be  modified  as 
follows: 


(1)  For  points  fixed  in  both  x  and  y  directions,  Ar  ;  Arv  -  C: 

x  y 


£  --  C  -  £  =  C  -  o 


(43a) 


(2)  For  points  fixed  in  the  x  direction  only,  Ar  s  0: 


r*  s 

Txvj 


f*  -  o 


(43b) 


(3)  For  points  fixed  in  y  direction  only,  =  0: 


r*  = 

txy 


C* 

+3* 


c  *  ^  - 


f  *  m  O 
i-33  V 

fijx  (*y 


(43c) 


For  the  structure  considered  in  this  report  all  boundary  nodal  points  were 
fully  fixed,  and  Eq.  43a  was  applied. 

Sequence  of  Operations 


Only  the  main  operations  of  the  computer  program  will  be  described; 
the  details  of  coding  will  be  omitted.  The  program  was  written  for  sn 
IBM  704  computer  with  32768  storage  locations.  The  operation  of  the  program 
is  flexible  in  that  both  input  and  output  can  be  "on-line"  or  may  be 
effected  "off-line"  through  the  use  of  magnetic  tapes  and  peripheral 
equipment . 

*nPut  Dnta  -  For  the  purpose  of  numerically  defining  a  structure,  all 
nodal  points  and  elements  are  numbered  as  illustrated  in  Fig.  7.  The 
numerical  description  is  read  into  the  machine  in  the  form  of  punched 
cards,  by  the  following  four  arrays: 

A.  Parameter  Array  (6  numbers) 

1.  Numbered  elements 

2.  Number  of  nodal  points 

3.  Number  of  boundary  points 

4.  Over-ralaxation  factor 

5.  Convergence  limit 

6.  Coefficient  of  thermal  expansion 

B.  Element  Array  (9  numbers  per  element) 

1.  Element  number 

2.  Number  of  nodal  point  i 

3.  Number  of  nodal  point  j 

4.  Number  of  nodal  point  k 

5.  Modulus  of  elasticity  E 

6.  Poisson’s  ratio  l) 


ASSUMED  DEFORMED  SHAPE 


FIG.  7  NUMBERING  SYSTEM  FOR  NODAL  POINTS  AND  PLATE  ELEMENTS 
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7.  Unit  weight  of  element '3c 

S.  Temperature  change  vithin„element  AT 
9.  Orthotropic  factor  =  ^ 

C.  Nodal  Point  Array  (7  numbers  per  nodal  point) 

1.  Nodal  point  number 

2.  x-ordinate 

3.  y-ordinate 

4.  X-load 

5.  Y-load 

6.  Initial  x-displacesent 

7.  Initial  7-displacement 

D.  Boundary  Condition  Array  (2  numbers  per  boundary  point) 

1.  Nodal  point  number 

2.  This  number  indicates  the  type  of  restraint;  "0”  for 
fixed  in  both  directions;  "1"  for  fixed  in  x-direction; 
and  "2”  for  fixed  in  y-direction. 

It  should  be  noted  that  for  a  fixed  boundary  point,  the  initial  displace¬ 
ment  is  the  final  displacement  of  the  point,  since  it  ±8  not  altered  by 
the  iteration  procedure. 

Formation  of  Element  Stiffness  Matrices  -  The  stiffness  matrix  for  each 
element  is  deteriri-.-i  f: —  ~q  *7  ?r  Eq  r>'~  The  basic  element 
dimensions  are  calculated  fro©  the  coordinates  of  the  connecting  nodal  points 


a j  -  xj  ~  xc 
bj  =  -  Vi 

h*  ~  y*  ~ 

where  i,  j,  and  k  are  the  nodal  point  numbers  of  the  three  connecting 
points  and  are  given  in  the  element  array. 

Formation  of  Complete  Stiffness  of  System  -  Because  of  the  large  matrices 
that  are  developed  in  a  plane  stress  system,  the.  formal  matrix  procedure 
represented  by  Eq.  1  is  not  used  directly.  Since  the  complete  stiffness 
matrix  contains  many  zero  elements,  only  the  non-zero  elements  are  developed 


3S 


ci  reliise:  by  the  pTOgrsr,  thee .  It  Is  possible  to  treat  Iirji  sisss 
stress  systeas  vitbset  egeeiit;  tbs  storage  capacity  -of  the  coapcter . 

Feme  tics  of  Lead  Matrix  -  Tie  lead  aazrix  i  2^  i*  ccep-csed  ef  lire  ioadi, 
dead  leads,  and  tesperatnre  leads-  The  ercs  ticca  wbicb  are  csed  te  deter- 
teire  these  leads  bare  beer  presetted  it  the  preriecs  sect let  of  this  report 
Fcmxclcc  of  Bfodal  Flexibilities  -  The  cecal  point  flexibilities  ere  deter- 
rsf. aed  frem  the  pruvlccsly  developed  stiffness  cc-sf fiejerts.  The  flexi¬ 
bilities  associated  with  the  eccrdxry  Dedal  points  are  codified  by  the 
appllcstlac  of  Eqs.  c3  as  received - 

Iterative  Sciatica  -  The  repeated  application  of  £q.  42  at  ail  cccal  rcin-tf 
cccatitates  the  iterative  prccedbre.  Tbs  roc  of  the  absolute  nagnitoda  c£ 
the  t-eba-.arred  fortes  at  ail  soeal  points  (giver  by  £<5.  41}  is  else  cco.- 
poted  for  each  cycle;  this  son.  when  con pared  to  the  convergence  licit, 
serves  as  a  cbedc  so  convergence  of  the  yrccodnre.  la  ell  analyses 
presorted  i£  this  report,  this  scst  vss  redoted  to  less  than  1/1  OSS  of  the 
value  obtained  it  the  first  cycle  of  iteraci.ee. 

Calculation  of  Fleaart  Stresses  -  Proa  the  nodal  point  dlsplaceaeots , 
sith  the  aid  cf  is.  23,  the  el  swot  stresses  ^  ,  5^,  f  r*4  ere 

calculated,  is  added  iafontsctoa  the  principal  stresses  3"  and  and 

directions  c?  are  also  calc^lated- 

Patent  Infomaticu  -  it  desired  polar s  in  the  iteration  procedure,  nodal 
displacements  arc  element  stresses  are  printed.  Fig.  3  illustrates  the 
£om  of  the  computer  output,  in  a  typical  case. 

Tinlcjt 

The  computational  tine  required  by  the  prog  ran.  is  approximately 
equal  to  0.07  seconds,  where  u  equals  the  number  01  nodal  points  and 
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c  equals  the  number  of  cycles  of  Iteration.  The  number  of  cycles  required 
depends  on  the  accuracy  Gf  the  initially  assumed  displacements  and  on  the 
desired  degree  of  convergence.  For  the  structure  considered  in  this  report, 
the  ccerputer  nice  per  solution  vas  approximately  7  minutes  for  the  coarse 
oeeh  and  17  minutes  for  the  fine  mesh  idealization.  Between  70  and  100 
cycles  of  iteration  were  used  in  typical  cases. 


SCHEDULE  0?  ANALYSES 


The  Structural  System 

The  dimensions  and  properties  of  the  structural  system  considered 
in  this  investigation  were  supplied  by  the  Contracting  Office,  Fig.  $ 
shows  the  basic  geometry  of  the  section  of  dam  which  was  analysed. 
Monolith  16.  This  was  one  of  the  cost  severely  cracked  portions  of  the 
structure.  The  extent  of  the  cracking  could  not  be  defined  exactly, 
but  it  appeared  that  there  was  one  principal  crack  which  extended  from 
near  the  center  of  the  base  to  roughly  two-thirds  of  the  height  of  the 
section,  as  shewn  in  the  figure.  The  physical  properties  assumed  for  the 
various  materials  in  the  system  are  listed  in  Table  I,  below. 


Table  I.  ASSUMED  PROPERTIES  OF  MATERIALS 


Concrete: 

Modulus  of  Elasticity 

Ec 

e 

2.0  x  106  psi 

Poisson* e  Ratio 

o 

- 

0.17 

Unit  Height 

4 

r 

150  pcf 

Them*!  Coefficient 

cX. 

— 

7.0  x  10  per 

Foundation  Rock: 

Modulus  of  Elasticity 

% 

z 

5.0  x  106  psi 

Poisson's  Ratio 

D 

= 

0.17 

Thermal  Coefficient 

aL 

z 

7.0  x  XO'6  per 

Vertical  Modulus 

% 

- 

1.0  x  106  psi 

(Orthotropic  cases) 

Water: 

Unit  Height 

£ 

z 

62.5  pcf 

Temperature  &T  *  '35°  in  body  of  dam  decreasing  to 

change: 

G°F  about  30  feet  below  the  surface  of  the 


foundation  rock  (See  Fig.  10} 


FIG.  9  BASIC  GEOMETRY  OF  MONOLITH 


The  relatively  low  modulus  of  elasticity  assumed  for  the.  concrete 
was  chosen  in  order  to  account  for  the  effects  of  creep  which  would  be 
expected  to  have  occurred  in  the  twenty  years  since  the  dam  was  built. 

Dead  weight  of  the  foundation  rock  was  not  considered  because  the  de¬ 
formations  and  stresses  resulting  from  this  effect  existed  prior  to  the 
construction  of  the  dam.  Only  conditions  resulting  from  the  construction 
of  the  dam  were  included  in  the  investigation. 

Cases  Considered 

The  principal  variables  considered  in  these  studies  included  the 
form  of  loading,  the  crack  height,  the  type  of  foundation  material 
(whether  isotropic  or  orthotropic)  and  the  possibility  of  uplift 
pressures  at  the  upstream  edge  of  the  base.  Several  factors  were 
considered  in  selecting  the  most  significant  combinations  of  these 
variables.  For  example,  temperature  changes  are  known  to  have  played 
an  important  part  in  the  behavior  of  the  structure,  and  are  assumed 
to  have  been  the  primary  cause  of  cracking.  Therefore,  temperature 
effects  were  included  with  each  of  the  basic  configurations.  However* 
because  the  ..  .  _1  s'.  a  ,:t  rep 

. .  it  was  decided  to  use  the  thermal 

analysis  primarily  for  determining  the  extent  of  the  crack  opening,  and 
to  calculate  separately  the  stresses  due  to  dead  and  live  loads.  Analyses 
of  an  uncracked  section  were  included  so  that  the  specific  influence  of 
the  crack  cn  the  stress  distributions  resulting  from  each  of  the  loading 
conditions  could  be  evaluated.  The  orthotropic  foundation  was  considered 
only  in  those  cases  believed  to  be  most  representative  of  actual  conditions 
in  order  to  determine  the  importance  of  this  factor.  Similar  considerations 
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dictated  the  choice  of  the  case  in  which  uplift  pressures  were  included. 
The  list  of  all  cases  studied  is  presented  in  Table  II, 


Table  II.  SCHEDULE  OF  ANALYSES 


Case 

grgck.H^ht. 

Foundation 

A 

None 

Dead  Load 

Isotropic 

B 

tt 

Dead  Live 

1? 

C 

it 

D  4*  L  4*  Temp. 

U 

D 

Two- thirds 

Dead  Load 

It 

E 

rt 

Dead  Live 

ft 

F 

it 

Dead  4*  Temp. 

rt 

G 

tt 

D  ^  L  4>  Temp. 

ft 

H 

tt 

Dead  4*  Live 

Orthotropic 

I 

Seven-Ninths 

Dead  Load 

Isotropic 

J 

tt 

Dead  4*  Live 

It 

K 

It 

Dead  4  Temp. 

:i 

L-l 

U 

Dead  4  Live 

it 

L-2 

ft 

ft  »» 

it 

M 

II 

n  tt 

Orthotropic 

H 

2/9  to  7/9 

»i  n 

Isotropic 

*0-1 

Two- thirds 

Temperature 

Isotropic 

*0-2 

tt 

Dead  4*  Live 

l> 

Uplift 

No 

ft 

5* 

tt 

rt 

it 

it 

ft 

ti 

»t 

ti 

3  nodal  points 
3%  '* 

No 

3  nodal  points 
No 

•  f 


*These  cases  were  includeu  to  provide  a  check  on  the  analyses  made  by 
Hr.  D.  McHenry  of  the  Portland  Cement  Association;  material  properties 
were  those  used  in  the  PCA  analyses  rather  than  those  cited  above. 
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Analysis  Procedure 

Foundation  Displacements  -  The  finite  element  approximations  used  in 
the  analysis  of  this  dam  are  shown  in  Fig.  10.  The  coarse  mesh  approxi¬ 
mation  was  employed  first  in  order  that  a  large  foundation  zone  might  be 
included  in  the  analysis.  With  this  idealization,  the  magnitude  of  the 
foundation  deformations  could  be  determined  with  good  accuracy.  Then 
the  displacements  calculated  along  the  dashed  line  in  the  coarse  mesh 
approximation  were  introduced  as  boundary  displacements  at  the  base  of 
the  fine  mesh  system.  By  this  procedure  it  was  possible  to  retain  the 
effect  of  the  complete  foundation  and  still  not  devote  a  large  number  of 
elements  to  its  treatment  in  the  final  analysis. 

Separate  coarse  mesh  analyses  were  carried  out  for  each  basic 
loading  condition:  dead  load,  dead  plus  live  load,  dead  plus  temperature, 
etc.  Each  set  of  calculated  boundary  displacements  was  then  used  in  the 
analysis  of  the  fine  mesh  case  for  the  corresponding  loading.  It  was 
assumed,  however,  that  the  extent  of  the  crack  would  not  affect  the 
boundary  displacements,  thus  the  same  displacements  were  introduced  for 
a  given  load  case  regardless  of  the  crack  condition.  It  should  be  noted 
that  water  loads  were  applied  to  the  appropriate  foundation  elements  as 
well  as  to  the  dam  structure  in  both  coarse  and  fine  mesh  live  load 
analyses . 

Treatment  of  Crack  *  The  presence  of  an  open  crack  introduces  no  compli¬ 
cation  in  the  finite  element  analysis  procedure;  each  side  of  the  crack 
is  treated  just  like  any  other  external  boundary.  Nodal  points  are 
defined  along  both  sides  of  the  crack  as  shown  in  Fig.  11,  and  displace¬ 
ments  of  the  two  sets  of  points  are  assumed  to  take  place  independently. 
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In  order  that  this  assumption  be  reasonable,  it  is  necessary,  of  course, 
that  each  pair  of  points  move  apart  rather  than  toward  each  other.  Thus 
it  was  assumed  that  the  dam  was  subjected  first  to  the  dead  load  plus 
temperature  change  conditions,  which  resulted  in  an  opening  of  the  crack. 
Then  if  the  subsequent  live  loading  did  not  completely  eliminate  the 
separation  at  any  pair  of  nodal  points  it  is  clear  that  the  two  sides 
would  continue  to  act  independently  as  assumed. 

Pplift  Pressures  -  It  was  recognised  at  the  outset  of  this  investigation 
that  there  might  be  a  tendency  for  tensile  stresses  to  develop  at  the 
upstream  edge  of  the  base  of  the  dam.  Because  the  tensile  strength  at 
this  point  is  essentially  nil,  it  would  be  expected  to  crack  open  under 
any  loading  condition  for  which  tensile  stresses  were  indicated.  Further¬ 
more,  since  this  crack  would  be  exposed  to  the  full  hydrostatic  pressure 
at  the  base  of  the  dam  (in  live  load  conditions),  uplift  pressures  would 
be  developed  over  the  full  extent  of  the  crack. 

The  uplift  pressure  analyses  were  treated  in  exactly  the  same 
manner  as  were  the  vertical  crack  analyses,  with  separate  nodal  points 
being  located  along  the  base  of  the  dam  and  in  the  foundation  rock,  as 
shown  in  Fig.  12.  The  only  difference  was  that  hydrostatic  live,  load 
forces  were  applied  at  the  nodal  points  of  this  crack,  acting  normal  to 
both  the  upper  and  lower  surfaces,  as  shown.  It  was  assumed  that  the 
crack  extended  downstream  a  distance  of  three  nodal  points  (about  32  feet) 
because  a  grout  curtain  and  relief  drain  system  located  in  this  vicioity 
would  prevent  the  uplift  pressure  from  extending  further.  (One  analysis 
was  made  in  which  the  uplift  pressure  wss  applied  over  a  distance  of 
about  42  feet,  but  this  condition  was  not  considered  realistic). 
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RESULTS  OF  ANALYSES 

The  stress  components  in  each  element  and  the  displacement  components 
of  each  nodal  point,  which  were  calculated  by  the  digital  computer  and 
printed  in  tabular  form  (see  Fig.  8),  constituted  the  direct  results  of 
the  analysis  procedure.  However,  the  mass  of  numerical  data  was  too  bulky 
to  be  digested  readily  in  this  fora,  and  a  considerable  effort  was 
expended  in  representing  the  most  significant  of  the  data  in  graphical 
plots. 

Stress  data  were  presented  in  two  different  forms,  stress  vectors 
and  stress  contours.  The  stress  vectors  are  merely  direct  graphical  plots 
of  the  princip'l  stresses  in  each  elemant.  They  ere  plotted  from  the 
center  of  the  element,  and  consist  of  arrows  showing  the  sense  and 
direction  of  the  stresses,  the  length  of  the  arrow  representing  the 
magnitude.  Study  of  the  stress  vector  diagrams,  which  are  compiled  in 
the  Appendix,  reveals  a  great  deal  of  detailed  information  about  the 
state  of  stress  in  the  structure.  However,  a  more  easily  interpreted 
representation  is  provided  by  the  stress  contours  which  are  presented 
and  discussed  in  the  first  section  of  this  chapter. 

Displacement  data  are  also  presented  in  two  different  sets  of 
figures.  The  boundary  displacements  are  merely  figures  depicting  the 
deformed  shape  of  the  cross-section  under  load.  These  are  of  qualitative 
interest,  but  have  no  great  significance  with  regard  to  the  performance 
of  the  structure.  The  crack  displacements  are  shown  in  these  figures, 

l 

but  in  order  to  provide  a  clearer  impression  of  the  extent  of  the  crack 
opening,  and  of  the  effect  on  it  of  the  various  load  combinations, 
separate  figures  have  been  constructed  showing  only  the  crack  displacements. 
The  displacement  figures  are  presented  and  discussed  in  the  latter  part 
of  this  chapter. 
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Stress  Contours 

Lines  connecting  points  of  equal  stress  intensity  in  the  cross-section 
are  called  iso-static  lines  or  stress  contours.  Such  contours  were  drawn 
in  each  case  for  the  vertical  normal  stress  (J^  ,  the  shear  stress  Tx^  , 

and  for  the  principal  stresses  CT ^  and  The  ®a1ior  principal  stress, 

0*1  is  defined  as  the  maximum  tensile  (or  minimum  compressive)  stress, 
while  the  minor  principal  stress,  (H,,  is  the  maximum  compressive  component. 
In  constructing  the  stress  contours,  it  was  assumed  that  the  stress  values 
computed  for  each  element  existed  at  the  centroid  of  the  element,  and  the 
contours  were  based  on  these  controls.  However,  a  certain  amount  of 
smoothing  was  required  in  converting  the  discontinuous  finite  element 
results  into  a  representation  of  the  actual  continuous  system. 

Also  presented  for  each  case  are  plots  of  the  0^  and  stress 

distributions  across  the  base  section  of  the  dam.  These  curvep  were 

defined  by  the  intersections  of  the  stress  contours  with  the  base  line. 

The  resultant  forces  represented  by  these  base  stress  distributions  were 

evaluated  directly  from  the  graphs  by  means  of  planimeters.  Comparisons 

of  the  calculated  force  resultants  with  the  applied  loads  for  all  cases 

are  presented  in  Table  III,  at  the  end  of  this  section. 

Case  A  (Fig,  13)--  The  base  stress  distribution  computed  In  this  uncracked 

section  is  very  nearly  of  the  simple  trapezoidal  form  assumed  by  the 

elementary  theory,  although  a  slight  tendency  toward  stress  concentration 

is  evident  at  the  heel  of  the  dam  (the  expected  effect  of  foundation 

elasticity).  The  shear  stress  distribution  across  the  base  for  this 

simple,  dead  load  condition  merely  represents  the  resistance  provided 

by  the  foundation  rock  to  the  Poisson* s  ratio  expansion  tendencies  in  the 

» 

concrete.  The  principal  stress  plots  show  that  no  tensile  stress  is 
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developed  at  any  point  in  the  section,  and  that  the  maximum  compressive 
stress  is  about  250  psi,  only  slightly  more  than  would  have  been  predicted 
by  elementary  theory. 

Case  B  (Fig.  14) —  Comparison  of  this  case  with  Case  A  shows  the  influence 
of  the  live  load  on  the  stresses.  Some  definite  shifts  in  the  distributions 
are  evident,  due  to  the  change  in  location  aud  direction  of  the  resultant 
applied  forces,  but  no  significant  change  has  occurred  in  the  maximum  stress 
values;  actually  the  maximum  stresses  have  been  reduced  slightly  by  the 
live  load.  The  most  notable  change  from  Case  A  is  the  reversal  of  the  base 
shear  stress  direction  in  the  region  of  the  heel.  This  is  a  direct  xwv.lt 
of  the  lateral  hydrostatic  force,  which  in  this  case  must  be  balanced 
by  the  base  shear  stresses. 

Case  C  (Fig.  15)—  In  this  figure  it  is  seen  that  the.  addition  of  the 
thermal  loads  causes  very  substantial  changes  in  the  stress  magnitudes  and 
directions,  especially  near  the  base  where  the  foundation  resist,  the  thermal 
strains.  The  primary  thermal  effect  is  a  lateral  contraction  of  the  con¬ 
crete  together  with  the  restraint  to  this  contraction  which  is  provided 
by  the  foundation  rock.  Although  the  assumed  thermal  gradient  in  the 
vertical  direction  tends  to  distribute  this  constraining  effect  over  a 
fairly  deep  portion  of  the  dam  near  the  base,  the  stresses  which  result 
are  quite  large. 

The  direct  consequence  of  the  base  restraint  is  a  system  of  tensile 
stresses  0^  in  the  horizontal  direction  (near  the  base)  accompanied  by 
shear  stresses  (  )  on  horizontal  planes.  These  shear  stresses  oppose 

those  due  to  dead  load,  but  when  combined  with  live  load  effects  they 
produce  a  maximum  shear  stress  of  nearly  400  psi  near  the  heel.  At  the 
same  time,  the  tensile  normal  stresses  (  (TJ  )  in  the  base  zone  tend  to 
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cause  a  rounding  of  the  base  plane  (convex  downward),  and  this  rounding 
tendency  is  resisted  by  vertical  normal  stresses  (0*3). 

The  resulting  tensile  stress  of  400  psi  (as  wall  as  the  major 
principal  stress  of  ever  500  psi)  in  the  vicinity  01  the  heel  demonstrates 
the  severity  of  the  loading  produced  by  the  temperature  changes.  Th;.~fi 
the  fact  that  the  structure  cracked  is  not  surprising. 

On  the  other  hand,  it  should  be  recognized  that  the  assumed  thermal 
stresses  probably  are  considerably  more  severe  than  those  which  actually 
occurred,  at  least  in  the  vicinity  of  the  faces  of  the  dam.  These  zones 
probably  never  achieved  the  high  temperatures  which  developed  in  the 
mid-section,  thus  the  thermal  stresses  near  the  faces  would  have  been  .. 
considerably  less  than  the  computed  values.  The  assumed  conditions- 
probably  are  quite  realistic  in  the  mid-section,  however,  and  the  computed 
tendency  for  cracking  in  this  area  should  be  representative  of  actual 
conditions. 

Case  D  (Fig.  16)--  The  comparison  of  this  case  with  Case  A  demonstrates 
the  relatively  minor  influence  which  the  crack  had  on  the  stresses  due 
to  dead  load  alone.  The  only  notable  effect  is  with  respect  to  the 
stresses  resulting  from  Poisson’s  ratio  expansion.  The  lateral  stresses 
due  to  this  effect  are  greatly  relieved  by  the  vertical  crack  through  the 
center  of  the  section,  thus  the  corresponding  shear  stresses  at  the  base 
of  the  section  are  also  reduced.  Minor  normal  stress  concentrations  appear 
at  both  upper  and  lower  enus  of  the  crack,  but  the  maximum  stresses 
associated  with  these  are  negligible. 

Case  E  (Fig.  17)--  The  influence  of  the  crack  on  the  stress  distribution 
due  to  live  load  is  much  greater  than  on  that  due  to  dead  load,  as  is 
indicated  by  a  comparison  of  Cases  E  and  B.  The  cracked  section  tends  to 
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resist  the  lateral  load  as  two  separate  systems,  connected  only  at  the  top. 
Thus  the  normal  stresses  (  Otj )  associated  with  overturning  are  much 
greater  in  the  upstream  block  than  they  were  for  the  solid  dam, reaching 
a  value  of  about  400  psi.  Similarly,  thj  shear  force  is  resisted  only  by 
the  upstream  block,  causing  a  significant  increase  in  shear  stress  as 
compared  with  the  uncracked  case. 

Sizeable  stress  concentrations  also  were  developed  at  the  upper  end 
of  the  crack  in  this  case,  due  to  the  action  of  the  zone  above  the  crack 
in  forcing  continuity  of  displacement  between  the  two  blocks.  Two 
distinct  actions  were  involved  in  maintaining  continuity.  The  lateral 
force  acting  on  the  upstream  block  pressed  it  against  the  restraint  of 
the  downstream  block,  causing  horizontal  compressive  stresses  (<J^  )  to 
develop  above  the  crack.  At  the  same  time,  the  bending  of  the  two  blocks 
tended  to  cause  the  contact  surface  of  the  left  block  to  move  downward 
while  the  corresponding  surface  of  the  right  block  tended  to  move  up. 

The  resulting  shearing  stresses  in  the  zone  above  the  crack  put  the 
upstream  face  of  the  crack  in  vertical  tension  (  0'3  )  and  the  downstream 
face  in  vertical  compression.  This  <3$  tensile  zone  in  the  upstream 
block  near  the  tip  of  the  crack  is  of  considerable  interest,  in  that 
horizontal  cracking  was  observed  at  this  location  in  the  prototype. 

The  corresponding  compressive  stress  in  the  downstream  block  is  consider¬ 
able  larger  in  magnitude,  of  course,  because  the  dead  weight  stress  in 
this  case  is  an  additive  rather  than  a  reducing  influence.  The  maximum 
compressive  stress  in  this  zone  is  no  larger  than  the  maximum  base  stress, 
however,  and  thus  is  not  of  particular  interest. 

Cage  F  (Fig-  18)—  Comparison  of  this  case  with  Case  D  shows  the  relative 
importance  of  the  thermal  loads  on  the  cracked  section.  The  thermal 
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effects  everywhere  are  similar  to  those  observed  for  the  uncracked 
section  (see  the  comparison  between  Cases  B  and  C)  when  due  allowance 
has  been  made  for  the  fact  that  the  crack  has  converted  the  system  into 
two  practically- independent  structures.  As  was  noted  in  the  discussion 
of  Case  C,  however,  the  assumed  thermal  conditions  are  not  entirely 
realistic,  and  a  better  impression  of  the  importance  of  the  crack 
i 8  provided  by  the  comparison  of  Cases  E  and  B,  in  which  only  live  and 
dead  loads  are  considered. 

Case  G  (Fig.  19)—  The  differences  observed  between  Cases  G  and  E  are 
equivalent  to  the  differences  between  Cases  F  and  D,  in  that  in  each 
case  the  differences  represent  the  effect  of  the  thermal  load  on  the 
cracked  section.  On  the  other  hand,  the  comparison  of  Cases  G  and  C 
demonstrates  the  influence  of  the  crack  on  the  stress  distributions 
due  to  identical  load  conditions.  However,  because  the  stresses  resulting 
from  the  temperature  changes  are  quite  large,  and  because  the  assumed 
thermal  conditions  are  a  rather  crude  approximation  of  the  actual 
conditions,  this  comparison  is  considered  to  be  less  significant  than 
that  discussed  above  for  live  and  dead  loads  only  (Cases  E  and  B). 

Case  H  (Fig.  20)—  The  influence  of  the  orthotropic  foundation  is  indicated 
by  a  comparison  of  Cases  H  and  E,  since  otherwise  these  two  cases  were 
identical.  J.t  is  apparent  from  the  similarity  of  the  stress  contours  for 
the  two  cases  that  the  vertical  foundation  flexibility  has  very  little 
effect  on  the  distribution  of  stress  in  the  upper  regions  of  the  structure. 
The  only  noticeable  effect  is  observed  near  the  base,  and  is  most  evident 
in  the  Otj  and  tKy  values  plotted  for  that  section.  The  reduced  vertical 
foundation  modulus  considered  in  this  case  makes  the  dam  behave  relatively 
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more  like  a  rigid  body.  Thus  the  computed  strec  distribution  is  more 
like  that  associated  with  a  rigid  dam  on  a  flexible  base,  in  which  stress 
concentrations  tend  to  develop  at  toe  and  heel.  However,  it  will  be 
noted  that  these  increased  stress  concentration  tendencies  in  Case  H  as 
compared  with  Case  E  ere  not  great;  in  fact  they  are  almost  within  the 
range  of  accuracy  of  the  method  of  analysis. 

Case  I  (Fig.  21)--  This  case  corresponds  exactly  with  Case  D,  the  only 
d i  fference  being  that  the  crack  has  been  extended  upward  for  one  more 
nodal  point  to  1/1  height.  Comparison  of  these  cases  demonstrates  that 
the  additional  crack  height  has  a  negligible  influence  on  all  aspects  of 
the  stress  distribution  due  to  dead  load  alone. 

Case  J  (Fig.  22)--  The  relationship  between  Cases  I  and  D,  which  was 
discussed  above,  applies  also  to  Cases  J  and  E,  except  that  in  these 
latter  cases  the  live  load  also  is  acting.  Although  the  effect  is  not 
large,  it  is  apparent  that  the  increased  crack  height  does  result  in 
a  definite  increase  in  certain  critical  live  load  stress  values.  Over 
most  of  the  cross-section,  the  similarities  between  the  two  cases  are 
quite  remarkable;  but  the  maximum  tensile  and  compressive  stresses  at 
the  base  show  a  significant  increase  in  the  7/9  crack  case,  as  does  the 
maximum  tensile  stress  at  the  upst'eam  face  of  the  crack.  This  finding 
is  not  surprising,  of  course;  since  these  critical  stress  values  are  a 
direct  consequence  of  the  crack,  it  is  to  be  expected  that  extending 
the  crack  would  result  in  an  increase  in  their  magnitude. 

Case  K  (Fig.  23)--  Another  indication  of  the  effect  of  extending  the  crack 
height  is  shown  by  the  comparison  of  cases  K  and  F.  As  would  be  expected 
in  the  absence  of  live  load,  the  difference  shown  by  this  comparison  (for 
dead  and  thermal  loads)  1b  negligible. 
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Case  L-l  and  L-2  (Fig.  24  and  25)--  The  effect  of  uplift  pressures  at 
the  heel  of  the  dam  is  shown  by  a  comparison  of  Cases  L  and  J.  The 
principal  effect  of  the  uplift  condition  is  the  upward  movement  of  the 
heel  of  the  dam,  and  the  consequent  rotation  of  this  upstream  portion 
of  the  dam.  This  rotation  causes  a  significant  increase  in  the  ecm- 
pressive  stress  at  the  base  of  the  dam  on  the  upstream  side  of  the  crack, 
and  also  causes  increases  in  the  stress  concentrations  near  the  upper 
end  of  the  crack. 

Cases  L-l  and  L-2  differ  only  in  the  extent  of  the  uplift  crack 
assumed  at  the  heel  of  the  section.  In  Case  L-l  it  extended  downstream 
for  three  nodal  points  (about  32  ft.),  in  Case  L-2  it  extended  3%  nodal 
points  (about  42  ft.).  The  additional  uplift  in  Case  L-2  caused  an 
increased  stress  condition  as  would  be  expected;  however,  Case  L-l  is 
considered  more  realistic  because  a  system  of  drains  and  a  grout 
curtain  provided  in  the  prototype  should  be  effective  in  limiting 
uplift  effects  to  about  the  first  three  nodal  points. 

The  maximum  stresses  observed  in  Case  L-l,  the  most  critical  case 
considered  (except  for  L-2),  amounted  to  over  400  psi  compression  at  the 
base  upstream  of  the  crack  and  on  the  downstream  face  of  the  upper  end 
of  the  crack.  At  the  same  time  a  tensile  stress  of  about  150  psi  was 
observed  on  the  upstream  face  near  the  upper  and  of  the  crack.  The 
tension  at  the  heel  of  the  dam  which  was  noted  in  Case  J  has  been  eliminated, 
of  course,  by  the  opening  of  the  crack  at  this  location,  and  in  its  place 
is  the  100  psi  hydrostatic  uplift  pressure.  Similarly,  the  shear  stress 
near  the  heel  has  been  eliminated  by  the  crack  and  the  entire  shear 
acting  on  the  upstream  block  is  resisted  by  the  remaining  uncracked 
portion  of  the  base. 
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Case  M  (Fig.  26)—  This  case  is  comparable  with  Case  H  in  that  it  represents 
the  condition  of  dead  plus  live  load  acting  on  the  system  with  the  ortho¬ 
tropic  foundation.  The  only  difference  between  this  and  Case  H  is  the 
extension  of  the  crack  to  the  7/9  height.  As  might  be  expected,  the  result 
of  this  extension  is  a  slight  Increase  of  the  stress  values  at  all  points 
of  concentration.  Of  greater  interest  is  the  comparison  with  Case  J, 
which  demonstrates  the  effect  of  the  orthotropic  foundation.  In  this 
comparison  it  is  quite  evident  that  the  softening  of  the  foundation 
tends  to  reduce  the  bending  constraint  acting  at  the  base  of  the  blocks, 
upstream  and  downstream  of  the  crack.  Thus  the  <1^  stress  gradient  across 
the  base  in  Case  M  is  considerably  less  than  in  Case  J.  This  softening 
affect  of  the  orthotropic  foundation  is  considerably  more  evident  here 
in  the  7/9  crack  height  case  than  it  was  with  the  2/3  crack  height 
(Case  H). 

Case  N  (Fig,  27)-“  In  this  case,  the  crack  has  been  assumed  to  extend 
upward  from  2/9  height  rather  than  from  the  base.  Otherwise  conditions 
are  identical  with  Case  L-l.  As  would  be  expected,  the  stress  concen¬ 
trations  here  are  greatly  reduced,  with  maximum  compressive  normal 
stresses  of  only  about  250  psi  rather  than  ^00  psi  as  in  Case  L-l.  In 
effect,  the  portion  of  the  dam  up  to  the  base  of  the  crack  acts  in  com¬ 
bination  with  the  foundation  as  a  sing?.e  unit.  The  renaining  portion 
which  is  subject  to  the  crack  effect  is  considerably  smaller  and  there¬ 
fore  is  stressed  less.  The  form  of  the  stress  distribution,  however,  is 
entirely  consistent  with  that  observed  before. 

Cases  0-1  and  0-2  (Figs.  2£.and  2 9)—  As  was  noted  previously,  these 
cases  were  included  to  provide  &  direct  comparison  between  results  obtained 
by  the  finite  element  method  and  those  computed  by  D.  McHenry  of  the 
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Portland  Cament  Association  using  the  lattice  anology.  The  properties 
assumed  by  McHenry  differed  significantly  from  those  employed  in  the 
analyses  described  above,  due  in  part  to  certain  limitations  in  his 
lattice  analogy.  Consequently  it  was  considered  desirable  to  employ 
these  same  property  assumptions  in  a  finite  element  analysis  for  the 
purpose  of  making  this  comparison.  Different  properties  were  assumed 
in  each  of  the  comparison  cases,  and  each  is  discussed  separately. 

Case  0-1  (Thermal  loading  only) : 

Temperature  change:  1  :  •  23°F  (uniform  in  dam,  no  change  in 

foundation) 

Moduli  of  Elasticity:  Foundation  E  =  5.0  x  10^  psi 

Dam  Ec=  1.90  x  106  psi 

Poisson's  Ratio:  V  »  0.33 

The  GjJ  stresses  computed  by  the  two  methods  are  seen  in  Fig.  28 
to  be  equivalent  for  all  practical  purposes.  The  stress  contours  are 
generally  the  same,  and  show  appreciable  differences  only  at  the  points 
of  stress  concentration  along  the  base  of  the  dam.  At  the  toe  and  the 
heel  of  the  section,  these  differences  are  due  primarily  to  the  boundary 
approximations  imposed  by  the  lattice  analogy  used  by  McHenry.  The 
structure  actually  analyzed  by  that  approach  is  one  in  which  the  sloping 
surfaces  of  the  dam  are  represented  by  steps  in  the  rectangular  mesh. 
Thus  at  each  step  a  stress  concentration  is  found  which  is  not  represen¬ 
tative  of  the  prototype  systan.  The  minor  differences  which  may  be 
seen  in  the  base  stress  concentrations  upstream  and  downstream  from  the 
crack  may  be  due  to  differences  in  the  methods  of  representing  the 
crack.  In  any  case,  the  substantial  agreement  between  results  obtained 
by  methods  as  different  as  these  is  good  evidence  of  the  reliability  of 


both  methods. 
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Case  0-2  (dead  and  live  load) : 

Modulus  of  Elasticity;  Dam  and  foundation  E  s  4.0  x  10^  psi 
Poisson’s  Ratio:  'O  -  0.33 

Unit  Weight:  Concrete  =  150  pcf 

The  principal  stress  distributions  for  this  loading  (Pig.  29)  are 
seen  to  be  remarkably  similar,  with  only  very  slight  differences  to  be 
noted  in  the  stress  concentrations  at  the  toe  and  the  heel.  These  are 
due  undoubtedly  to  the  step-wise  approximation  of  the  sloping  faces 
which  was  used  in  the  lattice  analogy,  and  are  of  no  practical  importance. 

Static  Checks  of  Force  Resultants 

In  applying  the  finite  element  method,  the  digital  computer 
automatically  satisfies  the  equilibrium  requirements  imposed  on  ^b.e 
structure;  in  other  words,  internal  forces  are  determined  which  exactly 
balance  the  external  loads  acting  at  each  nodal  point.  Thus,  there  is 
no  question  that  each  portion  of  the  finite  element  idealization  is  in 
equilibrium.  However,  a  further  stage  of  approximation  is  involved  in 
establishing  the  element  stresses,  and  particularly  in  constructing  the 
stress  contours  representing  these  element  stresses.  Thus  it  is  of  some 
interest  to  compare  the  resultant  forces  represented  by  the  calculated 
stress  distributions  with  the  external  load  resultants  applied  across  the 
same  sections.  Such  checks  were  made  at  tb"  base  of  the  dam  section  for 
all  cases  considered,  and  also  with  respect  to  the  small  triangular  block 
above  the  tip  of  the  crack  for  one  case. 

Base  Force  Resultants"  A  comparison  of  the  resultant  forces,  determined 

from  the  calculated  base  normal  and  shear  stress  distributions,  with 

the  total  vertical  and  horizontal  applied  load  components,  F  and  F 

y  3c 

respectively,  is  presented  in  Table  III. 
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TABL2  III.  '  CHECK  0?  BASE  FORCE  EQUILIBRIUM 


Case 

_S__ 

Error 

F, 

Error 

A 

3300X 

3350T 

-1.5% 

0 

0 

0 

B 

3740 

3650 

2.5 

1570 

1530 

2.6 

C 

3780 

3650 

3.6 

1640 

1530 

7-2, 

t 

D 

3400 

3350 

1.5 

0 

0 

0 

E 

3650 

3650 

0 

1640 

1530 

7.2 

F 

3250 

3350 

-3.0 

0 

0 

0 

G 

3800 

3650 

4.1 

1540 

1530 

0.7 

H 

3650 

3650 

0 

1620 

1530 

5.9 

I 

3350 

3350 

0 

0 

0 

0 

J 

3600 

3650 

-1.4 

1590 

1530 

3.9 

K 

3280 

3350 

-2.1 

0 

0 

0 

L-l 

3700 

3650 

1.4 

1520 

1530 

-0.6 

L-2 

3600 

3650 

"1.4 

1520 

1530 

-0.6 

M 

3700 

3650 

1.4 

1550 

1530 

1.3 

8 

3700 

3650 

1.4 

1560 

1530 

2.0 

In  general  it  is  believed  that  these  results  represent  a  very  satisfactory 
degree  of  agreement,  and  that  the  triangular  element  mesh  site  which  was 
employed  is  justified.  A  finer  definition  of  the  stresses  could  have 
been  achieved,  of  course,  by  using  smaller  elements,  but  at  a  considerable 
cost  in  computational  effort.  It  is  believed  that  no  greater  precision 
could  be  justified  for  the  purpose  of  this  project. 
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Block  Above  Crack"  From  the  complete  stress  distribution  information 
provided  by  these  analyses  it  is  possible  to  evaluate  the  equilibrium  of 
forces  acting  on  any  desired  portion  of  the  cross-section.  Of  particular 
interest  in  this  regard  is  the  small  triangular  block  above  the  tip  of  the 
crack,  which  transmits  all  of  the  forces  acting  between  the  two  basic 
portions  of  the  dam.  Fig.  30  shows  the  normal  and  shear  stress  distri¬ 
butions  in  this  zone  of  the  cross-section  for  Case  J,  the  case  with 
the  seven-ninths  crack  height  subjected  to  dead  plus  live  loads.  The 
normal  and  shear  stress  distributions  acting  on  the  vertical  and  horizontal 
surfaces  bounding  this  little  block  are  shown,  in  addition  to  their 
resultant  forces.  The  resultant  forces  acting  on  the  block  are  seen  to 
balance  to  within  about  4  per  cent,  which  is  a  very  &  at  ".s  factory  result, 
especially  considering  the  rather  small  number  of  elements  used  to  define 

the  state  of  stress  in  this  local  region. 

• 

Effect  of  Crack  on  Base  Stress  Distribution-  In  the  discussion  of  the 
stress  contours  above,  it  was  noted  that  the  effect  of  the  crack  on  the 
stress  distribution  was  indicated  directly  by  a  comparison  of  cracked 
section  results  with  the  uneraeked  section  results  for  like  loading 
conditions.  It  would  have  been  relatively  simple  to  arrange  the  program 
so  that  the  digital  computer  would  automatically  make  such  a  comparison, 
subtracting  stress  results  for  the  uncracked  cage  from  those  for  the  cracked 
section.  This  was  not  done  for  the  present  investigation  because  it  was 
considered  that  the  results  would  not  be  of  sufficient  interest  to  justify 
the  programming  effort.  However,  to  indicate  the  character  of  such  an 
analysis,  the  effect  of  the  crack  on  the  base  stress  distribution  for  a 
typieai  loading  was  evaluated  by  manually  taking  the  difference  between 
the  base  stresses  for  Cases  B  and  J.  The  results,  shown  in  Fig.  31,.  depict 
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the  base  normal  and  shear  stresses  caused  by  the  7/9  crack  when  the  section 
la  subjected  to  live  plus  dead  loading.  The  stress  concentrations  due  to 
the  crack,  amounting  to  a  maximum  of  about  200  psi  in  compression,  are  quite 
evident  in  this  plot. 

Boundary  Displacements 

The  displacement  vectors  of  each  of  the  boundary  nodal  points  have 
been  plotted  for  the  moat  significant  cases  which  were  treated.  Lines 
connecting  the  end  points  of  these  vectors  show  the  deformed  shape  of  the 
structure,  to  a  greatly  exaggerated  deformation  scale.  The  reference  base 
with  respect  to  which  all  displacements  have  been  computed  is  the  base  of 
the  coarse  mesh  foundation  elements  (Fig.  10),  thus  foundation  deformations 
as  well  as  those  of  the  concrete  structure  are  represented. 

Uncracked  Section  (Fig.  32  )-  The  only  notable  result  presented  in  this 
figure  is  the  relative  magnitude  of  the  thermal  deformations  as  compared 
with  those  due  to  dead  and  live  loads.  The  fact  that  serious  cracxing 
has  been  observed  in  the  structure  is  not  at  all  surprising  in  view  of 
these  very  large  thermal  strains. 

Two-Thirds  Crack  Height  (Fig.  33)-  Comparing  Figs.  32  and  33,  it  is  apparent 
that  the  only  effect  of  the  crack  on  the  boundary  deformations  is  a  very 
slight  increase  in  downstream  displacement  resulting  from  the  live  load. 

Thus  it  is  evident  that  the  cracked  structure  behaving  as  two  blocks  is 
slightly  more  flexible  than  the  uncrasked  system,  as  would  be  expected. 

The  introduction  of  the  orthotropic  foundation  leads  to  the  expected  large 
increase  in  the  vertical  displacement  of  the  structure „  but  has  no  other 
apparent  effect  on  the  boundary  displaeeaents. 
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Seven-Ninths  Crack  Height  (Fig.  34)-  Remarks  similar  to  the  above  apply 
equally  to  the  first  two  cases  presented  in  Fig.  34;  the  additional  crack 
height  considered  here  has  a  very  small  influence  on  the  boundary  displace¬ 
ment.  The  addition  of  uplift  pressures  at  the  heel  of  the  section  causes 
an  additional  downstream  rotation  of  the  dam,  as  is  shown  in  the  third 
sketch.  The  increased  stiffness  provided  by  closing  the  crack  for  the 
lower  two  nodal  points  is  demonstrated  by  the  reduced  downstream  dis¬ 
placement  to  be  noted  in  the  fourth  sketch  in  this  figure. 

Crack  Opening  (Fig.  35) 

The  displacements  of  the  upstream  and  downstream  faces  of  the  crack 
were  generally  of  more  interest  than  the  outer  boundary  iisplacements. 

In  the  first  place,  the  amount  of  opening  of  the  crack  was  used  as  a 
check  on  the  properties  assumed  for  the  concrete  (i.e.  its  thermal 
coefficient  and  modulus  of  elasticity)  as  well  as  on  the  assumed  tempera¬ 
ture  changes.  Because  the  opening  of  the  crack  rather  than  its  actual 
displacement  was  of  principal  interest,  only  the  horizontal  components 
were  plotted  in  this  figure,  relative  to  the  original  position  of  the 
crack. 

Two  sets  of  results  are  presented  in  this  figure,  those  for  the  2/3 
crack  height,  and  those  for  the  7/9  crack  height.  For  each  system,  the 
crack  opening  due  to  temperature  alone  was  determined  by  taking  the 
difference  between  dead  load  plus  temperature  and  the  dead  load  (alone) 
results.  The  maximum  width  of  the  opening  due  to  temperature  was  found 
to  be  over  1/8  inch  for  the  2/3  crack  height  case  and  over  3/16  of  an  Inch 
for  the  7/9  case,  the  increase  in  opening  due  to  the  vertical  extension 
of  the  crack  being  quite  notable.  These  values  decreased  to  about  1/8 
inch  and  just  under  3/16  inch  respectively  when  dead  load  was  included, 
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a  more  realistic  condition,  of  course.  Addition  of  the  live  load  then 
caused  both  a  large  downstream  displacement  of  the  crack,  as  well  as 
tendencies  toward  closing  near  the  upper  end  (where  the  base  constraint 
was  less  effective).  The  amount  of  closure  was  somewhat  greater  in  the 
7/9  case,  but  even  here  it  did  not  completely  close  at  any  point. 

However,  when  the  uplift  pressures  were  added,  it  was  found  that 
the  crack  did  actually  close  at  the  top  nodal  point  in  the  7/9  case. 

The  calculations  show  an  excess  of  closure  of  0.004  inches,  thus  the 
displacements  actually  would  have  been  slightly  less  than  these  computed 
values,  and  results  obtained  in  this  case  may  be  considered  slightly 
conservative.  In  other  words,  contact  pressures  across  the  upper  nodal 
point  of  the  crack  would  have  carried  some  of  the  stress,  and  results 
would  revert  back  slightly  toward  those  obtained  for  the  2/3  crack 
height  case. 

It  should  be  noted,  however,  that  only  in  the  most  severe  conditions 
assumed  did  any  crack  closing  develop,  and  since  the  crack  probably 
does  not  extend  all  the  way  from  the  base  to  7/9  height,  it  is  probable 
that  no  closure  actually  has  been  produced  on  the  prototype  structure. 


78 


C05CLUSI0HS 

The  results  obtained  in  this  investigation  permit  several  conclusions 
to  be  drawn  with  regard  to  a  number  of  different  aspects  of  the  study. 

It  should  be  pointed  out  that  the  investigators  were  not  requested  to 
establish  the  adequacy  of  the  structure  and  no  conclusions  of  this  type 
will  be  reported  (though  it  should  also  be  mentioned  that  none  of  the 
results  indicate  an  overstressed  or  unsafe  condition).  Specific  conclusions 
are  discussed  consecutively,  below. 

Method  of  Analysis 

One  of  the  most  important  results  derived  from  the  study  was  the 
damonstration  of  the  effectiveness  and  versatility  of  the  finite  element 
method  and  of  the  digital  computer  program  which  was  developed  to 
implement  the  method.  Stress  distributions  and  deformations  were  obtained 
for  a  variety  of  complex  plane  stress  systems  -subjected  to  some  rather 
awkward  loading  conditions.  The  input  data  required  to  define  any 
problem  is  easily  prepared;  the  most  time  consuming  part  of  the  task  is 
the  choice  of  the  most  suitable  finite  element  mesh  and  the  evaluation 
of  the  coordinates  of  the  nodal  points  of  the  system.  When  the  input  data 
has  been  prepared  properly,  the  operation  of  the  program  is  automatic 
and  the  analysis  can  be  carried  to  any  suitable  degree  of  convergence. 
Computer  time  requirements  are  quite  reasonable,  even  for  the  treatment 
of  complex  cases  such  as  these. 

The  comparison  of  this  method  with  the  lattice  analogy  employed  by 
D.  McHenry  demonstrated  that  both  methods  give  reliable  results.  The 
lattice  analogy  idealization  utilized  a  much  finer  mesh  than  was  considered 
desirable  in  the  finite  element  system.  In  any  analysis  of  this  type  it 
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is  necessary  Co  weigh  the  advantage  of  increased  precision  of  the  results 
against  the  cost  in  computer  time,  which  increases  rapidly  with  the 
number  of  elements.  In  general,  it  is  best  to  use  only  enough  elements 
to  obtain  the  necessary  accuracy;  this  is  particularly  true  when  a 
large  number  of  cases  are  to  be  treated.  Comparisons  of  the  computer 
time  requirements  of  the  two  methods  are  relatively  meaningless  because 
the  IBM  704  used  in  this  study  is  a  much  faster  and  more  powerful  machine 
than  the  LGP30  used  by  McHenry,  though  it  may  be  noted  that  the  analysis 
of  a  single  case  by  McHenry  took  about  200  times  as  long  as  to  do  a 
similar  problem  on  the  IBM  704. 

The  most  obvious  advantage  of  the  finite  element  method  over  the 
lattice  analogy  is  the  ease  with  which  it  can  account  for  arbitrary 
boundary  conditions (e.g.  sloping  surfaces),  arbitrary  material  properties 
(Poissons*  ratio  must  be  taken  as  1/3  in  the  lattice  analogy)  and 
arbitrary  thermal  variations.  Such  special  conditions  as  the  orthotropic 
foundation,  which  would  be  quite  difficult  to  represent  in  the  lattice 
analogy,  also  may  be  included  without  difficulty. 

Effect  of  Crack 

The  principal  reason  for  initiating  this  investigation  was  tho  question 
as  to  how  a  centrally  located  crack  would  affect  the  stress  behavior  of  a 
gravity  dam.  These  studies  showed  clearly  that,  although  several  rather 
significant  stress  concentrations  might  be  associated  with  the  crack, 
the  magnitudes  of  the  total  stresses  were  not  excessive  for  the  section 
geometry  considered  here,  being  less  than  500  psi  in  all  cases.  The 
principal  stress  concentration  zones  were  at  the  upper  and  lower  ends 
of  the  crack  and  at  the  heel  of  the  dam.  Tensile  stress  was  developed  at 
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the  heel  of  the  cracked  section,  and  at  the  upstream  face  of  the  crack 
near  its  upper  end.  It  is  of  interest  to  note  that  cracking  was  observed 
in  the  prototype  at  this  latter  location,  and  it  is  presumed  to  be 
present  in  the  vicinity  of  the  heel  also. 

Thermal  Loading 

The  finite  element  analysis  described  in  this  report  indicates  stresses 
due  to  thermal  changes  which  are  consistent  with  the  assumed  behavior  of 
the  structure.  However,  it  is  believed  that  the  magnitude  of  the  com¬ 
puted  stresses,  which  exceeds  500  psi  near  the  heel  of  the  cracked  section 
or  350  psi  at  the  base  of  the  section  with  the  crack  extending  to  7/9  height, 
is  considerably  greater  than  the  actual  thermal  stresses  developed  in 
these  initial  zones  because  several  of  the  assumptions  on  which  the  analysis 
is  based  are  not  satisfied. 

The  principal  differences  between  the  assumed  and  the  actual  conditions 
are  as  follows: 

1)  It  was  assumed  that  there  was  no  horizontal  temperature  gicdient, 
whereas  it  is  known  that  the  temperature  rise  (and  subsequent 
cooling)  must  be  significantly  smaller  near  the  exposed  surface 
of  the  section.  Tensile  stresses  near  the  surface  therefore 
would  be  significantly  less  than  computed, 

2)  It  was  assumed  that  the  entire  section  was  subjected  to  dead 
weight  and  cooling  effects  as  a  unit,  whereas  the  structure 
actually  was  built  over  an  extended  period  of  time,  and  cooling 
deformations  of  the  earliest  construction  took  place  before  the 
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final  sections  were  built.  Resistance  to  thermal  deformation 
(thermal  stress)  must  therefore  be  less  than  computed. 

3)  It  was  assumed  that  the  section  acted  as  a  homogeneous,  isotropic 
elastic  system  except  for  the  specified  vertical  and  uplift 
cracks,  whereas  it  is  known  that  many  additional  smaller  cracks 
developed  in  the  structure,  particularly  near  the  base.  These 
additional  cracks  would  have  the  effect  of  relieving  tensile 
stresses  wherever  the  computed  tensile  stresses  were  high 
enough  to  induce  cracking.  Specifically,  it  would  be  expected 
that  the  horizontal  (0^)  tensile  stress  over  the  entire  base 
section,  and  the  vertical  (CTy)  tensile  stress  near  the  heel 
would  be  reduced  by  this  effect. 

Because  of  these  several  factors  tending  to  reduce  thermal  stess, 
it  is  believed  that  the  computed  thermal  stresses  are  significantly 
higher  than  the  actual  values  in  the  critical  stress  regions.  Therefore 
it  was  deemed  advisable  to  base  the  principal  stress  analysis  conclusions 
of  the  study  on  dead  and  live  load  effects,  and  to  take  account  of  the 
thermal  analysis  with  regard  to  its  deformation  effects  only.  This  is  • 
equivalent  to  assuming  that  all  thermal  stress  was  dissipated  by  creep 
and  cracking. 

Orthotropic  Foundation 

Although  the  orthotropic  foundation  was  found  to  contribute  to  large 
vertical  displacements  of  the  section,  the  effect  of  the  orthotropy 
on  the  stress  distribution  was  Inconsequential.  Thus  any  question  as  to 
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whether  the  assumed  foundation  properties  are  correct  or  realistic  is 
of  littli*.  importance;  it  seems  clear  that  the  results  don't  depend  upon 
this  factor. 
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SUPPLEMENTARY  INVESTIGATION 

INTRODUCTION 

Objectives  of  the  Supplementary  Investigation 

During  the  course  of  the  original  investigation  regarding 
the  stress  distribution  in  Norfork  Dam,  a  number  of  topics 
came  up  which  were  of  interest,  but  were  beyond  the  scope  of 
the  contract.  Accordingly,  the  contract  was  modified  to  permit 
the  extension  of  the  investigation  into  two  of  these  areas, 
and  this  supplementary  report  describes  the  work  done  under 
the  extended  contract.  This  supplement  is  intended  merely 
as  an  addition  to  the  original  report  and  is  not  complete  in 
itself. 

One  of  the  topics  considered  here  involves  a  simple 
extension  of  the  analyses  performed  in  the  original  work,  to 
include  a  case  not  previously  treated.  In  this  new  case, 
which  is  designated  N~2,  the  central  crack  is  assumed  to  ex¬ 
tend  from  the  first  to  the  seventh  nodal  points  (i.e.,  from 
1/9  to  7/9  height).  Thus  it  lies  midway  between  cases  L— 1 
and  N  of  the  original  investigation.  It  is  now  believed  that 
this  new  case  is  closer  to  the  actual  conditions  in  the  dam 
than  any  of  the  previous  cases.  Case  N-3  is  identical  with 
N-2  except  that  the  reservoir  level  has  been  set  at  573  ft. 
elev.  rather  than  584  ft  as  was  assumed  for  all  other  live 
load  cases.  It  represents  the  state  of  stress  developed  by 
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the  maximum  water  level  yet  permitted  in  the  reservoir. 

The  other  subject  of  this  supplementary  study  is  the  case 
of  a  crack  extending  through  to  the  upper  surface  of  the  dam. 
Examination  of  the  dam  indicates  that  no  such  crack  exists; 
however,  this  condition  is  of  interest  as  a  limiting  case  as 
well  as  a  subject  of  scientific  curiosity.  The  physical 
properties  of  the  dam  were  assumed  the  same  for  this  case  as 
before,  and  a  similar  finite  element  idealization  was  adopted. 

As  before,  the  thermal  changes  were  used  merely  tc  establish  the 
initial  crack  opening  (under  dead  load  plus  thermal  conditions). 
The  thermal  stresses  were  considered  to  be  unrealistic,  due 
to  creep  effects,  and  were  not  included  in  the  final  stress 
distributions.  Isotropic  foundation  conditions  were  assumed 
in  all  of  the  new  analyses,  and  the  calculations  were  carried 
out  both  with  and  without  uplift  pressures. 

The  principal  difference  between  these  new  cases  and  the 
previous  analyses  is  that,  with  the  crack  extending  to  the 
surface,  the  section  becomes  so  flexible  that  addition  of  the 
live  load  causes  the  upstream  block  to  deflect  into  contact 
with  the  downstream  block.  This  "crack  closing"  condition 
makes  the  structure  non-linear  and  requires  a  major  modi¬ 
fication  of  the  analysis  procedure  in  order  to  account  for  it. 
Thus,  an  additional  result  of  this  supplementary  study  was  the 
demonstration  of  a  procedure  for  the  analysis  of  non-linear 
problems  by  the  finite  element  method. 
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Cases  Considered 

A  complete  listing  of  the  analyses  considered  in  this 
supplementary  investigation  is  presented  in  Table  Ha,  below. 

Table  Ila  SUPPLEMENTARY  SCHEDULE  OF  ANALYSIS 
Case  Crack  Heights _ Loading _ Uplift  Fig.  No. 


N-2 

1/9  to 

7/9 

Dead 

+ 

Live 

32  ft. 

36 

N-3 

1/9  to 

7/9 

Dead 

4- 

Live  (573ft.) 

32  ft. 

3  6a 

P-1 

Full 

Dead 

+ 

Temp 

None 

37  &  38 

P-2 

Full 

Dead 

+ 

Live 

None 

39 

P-3 

Full 

Dead 

+ 

Temp 

32  ft. 

P-4 

Full 

Dead 

+ 

L?  ve 

32  ft. 

40 

Q-l 

1/9  to 

full 

Dead 

4- 

Temp 

None 

Q-2 

1/9  to 

full 

Dead 

4- 

Live 

None 

41 

Q-3 

1/9  to 

full 

Dead 

+ 

Temp 

32  ft. 

Q-4 

1/9  to 

full 

Dead 

■t- 

Live 

32  ft. 

42 

Cases  P  &  Q  involve  the  new  full  crack  height  configura¬ 
tion  in  which  the  crack  closed  under  live  loading.  The  1st 
and  3rd  analysis  for  each  case  considered  dead  plus  thermal 
stresses  in  order  to  establish  the  initial  crack  opening. 

The  stresses  due  to  thermal  changes  were  not  considered 
significant,  however,  as  was  noted  above,  and  stress  contours 
are  presented  only  for  the  2nd  and  4th  analysis  for  each  case. 

METHOD  OF  ANALYSIS  FOR  CRACK  CLOSING 
General  Comments 


If  a  gravity  dam  section  is  cracked  all  the  way  to  the 
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surface,  the  crack  will  open  over  its  entire  length  under  the 
action  of  gravity  loads  and  thermal  changes.  As  the  reservoir 
is  filled,  however,  the  increasing  hydrostatic  pressures  will 
gradually  force  the  crack  closed  at  the  top;  and  as  more  water 
is  added  the  crack  will  close  further  and  further  down.  As 
the  crack  closes  and  contact  is  made  between  the  two  sides, 
the  behaviour  of  the  dam  changes.  Opposite  sides  of  the 
crack  are  no  longer  free  to  move  independently;  in  fact,  because 
of  the  roughness  of  the  crack  surface,  it  may  be  assumed  that 
the  two  sides  lock  together  upon  contact  and  subsequently 
displace  as  a  continuous  medium. 

This  changing  character  of  the  structure  with  crack 
closing  makes  it  infinitely  non-linear  in  reality.  However, 
in  the  finite  element  idealization  only  a  finite  number  of 
nodal  points  can  ccme  into  contact  along  the  crack  surface. 

Thus  an  incremental  approach  to  the  non-linear  problem 
becomes  feasible.  The  reservoir  is  filled  by  increments  and 
each  time  a  nodal  point  closes,  they  are  locked  together  and 
the  modified  structure  is  subjected  to  further  increments 
of  load.  For  the  purpose  of  this  analysis,  a  pair  of  nodal 
points  facing  each  ether  across  the  crack  (see  Fig.  11,  Final 
Report)  is  assumed  to  have  closed  when  their  relative  horizontal 
displacement  under  live  load  just  equals  their  initial  hori¬ 
zontal  displacement  due  to  dead  load  plus  thermal  changes. 
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The  stresses  and  displacements  developed  in  each  element 
of  the  dam  due  to  each  increment  of  loading  may  be  super¬ 
imposed  to  obtain  the  total  stress  in  the  non-linear  system 
at  any  stage  of  loading.  Thus  the  treatment  cf  the  crack 
closing  problem  reduces  to  the  analysis  of  a  succession  of 
varying  structures  due  to  load  increments  and  t"he  super¬ 
position  of  the  stresses  (and  displacements;  calculated  for 
each  increment. 

Uplift  pressures,  inhere  considered  in  these  analyses,  were 
assumed  to  act  ever  the  first  3  nodal  points  cf  the  base  ^32 
feet)  and  were  applied  at  the  full  hydrostatic  pressure 
appropriate  to  each  load  increment.  It  should  be  noted  that 
the  extent  cf  opening  of  the  uplift  pressure  crack  was  no_t 
assumed  to  increase  incrementally,  it  was  active  over  its 
full  extent  throughout  the  leading. 

Because  of  the  incremental  nature  ox  the  loading  applied 
to  +  he  dain  in  these  studies,  it  was  not  feasible  to  establish 
initial  foundation  displacements  with  a  coarse  mesh  system  as 
has  been  done  previously.  For  this  reason  a  larger  foundation 
zone  was  included  with  the  dam  section  in  order  that  the 
foundation  deformations  might  be  treated  realistically. 

As  was  noted  earlier,  the  theima!  stresses  computed  in 
establishing  the  initial  crack  displacement  condition  of  the 
structure  were  not  considered  to  be  realistic.  Thus  these 
stresses  were  stored  during  the  incremental  live  load  analysis 
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and  then  subtracted  off  prior  to  printing  the  final  stress  dis¬ 
tribution.  This  is  equivalent  to  assuming  that  the  thermal 
stresses  were  eliminated  by  creep  and  plastic  flow. 

The  Computer  Program 

The  plane  stress  program  used  in  all  the  previous  analyses 
formed  the  basis  of  the  crack-closing  analysis  program.  How¬ 
ever,  it  was  modified  by  the  addition  of  sub-routines  which 
accomplished  the  incremental  loading  and  structural  modifica¬ 
tion  operations  described  above.  The  basic  flow  sequence  of 
the  modified  analysis  is  as  follows 

1)  The  dead  plus  thermal  loading  case  is  solved,  for  the 
section  cracked  to  the  surface. 

2)  The  displacements  and  stresses  due  to  this  loading 
are  stored  and  the  loading  is  erased. 

3)  After  checking  that  the  water  level  is  not  at  the 
prescribed  maximum,  the  loads  due  to  a  specified 
increment  of  water  level  are  computed. 

4)  The  structure,  in  its  existing  crack  configuration,  is 
analyzed  for  the  incremental  water  leads. 

5)  The  resulting  horizontal  crack  displacements  are  checked 
against  the  previously  existing  crack  openings  to 

see  if  closure  develops  at  any  nodal  point. 

6a)  If  no  closure  occurs,  the  stresses  and  displacements 
are  superposed  on  the  previously  existing  results,  the 
loads  are  erased,  and  the  procedure  is  repeated  from 
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Step  3. 

6b)  If  exact  closure  occurs  at  any  nodal  point,  the  struct¬ 
ure  is  modified  by  altering  element  connections  so 
that  the  structural  form  is  that  of  a  section  with 
the  crack  height  lowered.  Then  stresses  and  displace¬ 
ments  are  superposed  into  storage,  the  loads  are 
erased,  and  the  procedure  is  repeated  from  Step  3. 

6c)  If  overdosing  occurs  at  any  nodal  point,  the  loads 
are  scaled  so  as  to  cause  exact  closing*  and  the 
structure  is  modified  as  above.  The  remaining  pro¬ 
portion  of  the  loads  is  saved  for  superposing  -with  the 
n en  load  increment.  The  scaled  stresses  and  displace¬ 
ments  are  superposed  into  storage,  and  the  procedure 
is  repeated  from  Stop  3. 

7)  When  the  maximum  water  level  is  reached,  the  prevailing 
stresses  and  displacements  in  storage  are  printed. 

8)  Finally  the  stresses  due  to  thermal  changes  are  sub¬ 
tracted  and  the  residual  stresses  (due  to  dead  plus 
live  load)  are  printed. 

Because  essentially  the  same  basic  program  was  used  for 
each  increment  of  loading,  the  analysis  of  ^  single  case 


*It  will  be  noted  that  this  scaling  of  loads  is  not  exactly 
equivalent  to  taking  a  smaller  water  level  increment.  However, 
for  small  increments  of  water  level,  the  difference  is  negligible. 
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involving  crack  closing  was  equivalent  to  approximately  10 
analyses  of  the  previous  type.  However,  the  high  speed  IBM 
7090  digital  computer  became  available  for  these  analyses, 
and  in  addition  certain  modifications  were  made  in  the  basic 
plane  stress  program  which  permitted  faster  operation.  As  a 
result,  a  complete  crack  closing  case  -  involving  approximately 
1200  cycles  of  iterations  in  all  -  could  be  solved  in  about  12 
minutes. 

RESULTS  AND  DISCUSSION 

Cases  N-2  and  N-3 

As  was  noted  earlier  these  cases  are  intermediate  between 
Cases  L-l  and  N  of  the  original  investigation  report,  in  that 
the  crack  extends  from  1/9  height  to  7/9  height  rather  than 
from  the  base  and  from  2/S  height  respectively  in  the  previous 
studies.  The  stress  contours  for  Case  N-2,  presented  in  Fig.  36, 
show  that  the  maximum  stresses  fall  between  the  values  given 
for  the  previous  studies  in  Fig.  24  and  27  respectively.  Thus 
the  compressive  stresses  at  the  upstream  side  of  the  base  of 
the  crack  are  reduced  from  L-l  but  are  greater  than  N,  while 
similar  effects  are  observed  with  respect  to  the  tensile 
stresses  near  the  top  of  the  crack. 

The  stress  contours  for  Case  ^-3  (shown  in  Fig.  36a)  are 
so  similar  to  Fig.  36  that  differences  can  be  noted  only  by 
careful  scrutiny.  The  loading  for  this  case  is  reduced  by  about 
10  per  cent  from  Case  N-2,  and  the  apparent  changes  in  stress 


93 


are  of  this  order  of  magnitude. 

Crack  Closing  Cases  (P  and  Q) 

Because  the  stress  distribution  of  the  systems  subjected 
to  crack  closing  depends  so  directly  upon  the  extent  of  crack 
closing,  ’t  is  advantageous  to  consider  just  the  crack  dis¬ 
placement  history  during  the  incremental  loading.  These  results 
are  presented  in  Fig.  37  for  Case  P  (cracked  to  the  base)  and 
in  Fig.  38  for  Case  Q  (cracked  to  1/9  above  the  base). 

Case  P-2,  in  which  no  uplift  is  considered,  shows  a 
regular  relative  displacement  effect  of  the  two  sides  of  the 
crack  as  the  load  is  increased,  with  a  final  condition  in  which 
the  crack  is  closed  from  the  top  to  7/9  height.  Thus  its 
final  configuration  is  quite  similar  to  Case  J  (Fig.  22). 

Case  P-4,  which  includes  the  uplift  crack  is  considerably 
more  flexible  (due  to  the  redirection  of  base  constraint) .  In 
fact,  in  the  initial  dead  load  plus  thermal  change  condition, 
the  analysis  indicates  that  the  uplift  crack  would  develop 
negative  displacements,  which  accounts  for  the  wide  opening  of 
the  vertical  crack  for  this  condition.  This  is  an  impossible 
situation,  of  course,  but  as  Jive  load  is  added  the  upstream 
block  tilts  rapidly  and  long  before  the  top  of  the  vertical 
crack  closes,  the  uplift  crack  has  opened.  Thus  the  analysis 
is  correct  for  an  initial  condition  consisting  of  dead  plus 
thermal  plus  partial  water  loading,  and  is  valid  for  all  sub¬ 
sequent  stages  of  loading.  The  greater  flexibility  of  the  P-4 
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case  (as  compared  *ith  P-2)  is  indicated  by  tne  fact  that  closure 
developed  from  the  top  to  5/9  height  of  the  crack. 

The  principal  difference  in  crack  displacements  shoun  for 
Case  0  in  Fig.  38  is  due  to  the  increased  stiffness  caused  by  the 
base  of  the  crack  starting  at  1/9  height.  This  additional 
connection  at  the  base  of  the  crack  inhibits  its  initial  opening 
under  dead  plus  thermal  loading.  Subsequent  live  loading  then 
is  able  to  cause  greater  closure  than  developed  in  Case  P-2,  i.e. 
from  the  top  to  6/9  height.  The  uplift  crack  produced  approxi¬ 
mately  the  same  result  in  Case  Q-4  as  in  Case  P-4,  including 
both  the  overdosing  effect  of  the  uplift  crack  before  live 
loading  (to  a  reduced  degree)  and  the  closing  of  the  vertical 
crack  from  the  top  to  5/9  height  under  full  live  load. 

The  stress  contours  for  Cases  P-2,  P-4,  Q-2,and  Q-4  are 
shoun  in  Figs.  39,  40,  41  and  42  respectively.  Comparison  of 
the  results  for  Case  P-2  dth  Case  J  (Fig.  22)  is  of  consider¬ 
able  interest  in  that  the  final  crack  configuration  is  the  same 
in  both  cases  but  vertical  adjustment  took  place  in  Case  P-2 
before  the  crack  closed  at  the  top.  Consequently  the  shearing 
forces  transmitted  across  the  crack  from  7/9  height  to  the  top 
are  much  less  in  this  case  than  for  the  corresponding  zone  in 
Case  J.  As  a  result,  the  live  load  is  carried  more  by  the 
upstream  block  in  Case  P-2  with  corresponding  increases  in 
the  stress  concentrations  at  the  base  of  the  block  (exceeding 
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600  psi  compression  and  300  psi  tension). 

Addition  of  the  uplift  pressures  in  Case  P-4  eliminates 
the  tensile  stress  concentrations  at  the  heel  of  the  section, 
of  course,  and  permits  a  greater  part  of  the  live  load  to 
be  transferred  to  the  downstream  block.  The  compressive 
stress  concentration  at  the  base  is  seen  to  be  reduced  from 
Case  P-2  (to  less  than  600  psi). 

Similar  results  are  evident  in  Cases  Q-2  and  Q-4,  except 
that  having  the  crack  closed  to  the  first  nodal  point  causes 
a  much  more  favorable  distribution  of  load  between  the  two 
segments  of  the  structure.  In  Case  Q-2,  the  maximum  tensile 
and  compressive  stresses  at  the  critical  horizontal  section 
are  only  about  150  and  450  psi  respectively.  Addition  of  the 
uplift  forces  in  Case  Q-4  again  reduces  the  maximum  compressive 
stress  at  the  critical  section  to  only  about  350  psi. 

One  of  the  basic  assumptions  made  in  the  "crack  closing" 
analyses  was  that  the  two  sides  of  the  crack  locked  together 
upon  contact,  and  that  no  relative  sliding  of  the  two  surfaces 
was  permitted.  This  assumption  is  valid,  of  course,  only  if 
the  coefficient  of  friction  between  the  surfaces  is  sufficient 
to  prevent,  slip.  In  order  to  assess  these  frictional  force 
requirements  the  shear  and  normal  stresses  developed  along 
the  part  of  the  crack  which  had  closed  were  evaluated  for  Cases 
P-4  and  Q-4.  The  results,  shown  in  Fig.  43,  demonstrate  that 
the  normal  and  shear  stresses  average  about  50  psi  in  Case  P-4 
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and  about  30  psi  in  Case  Q-4*  In  both  cases  the  shear  and 
normal  stresses  are  nearly  of  equal  value  thus  the  required 
coefficient  of  friction  is  about  1.  For  the  low  stress  values 
and  rough  concrete  surface  involved  here  this  seems  quite 
reasonable  and  the  assumed  behavior  is  believed  justified* 

CONCLUSIONS 

The  principal  contribution  of  this  supplementary  investi¬ 
gation  was  the  development  of  the  incremental  analysis  procedure 
which  made  possible  the  treatment  of  crack  closing  conditions* 
Although  the  incremental  loading  increases  the  computational 
effort  by  a  factor  of  about  10*  with  the  present  high  speed 
program  and  digital  computer  such  problems  may  be  handled  easily. 
Data  preparation  for  an  incremental  analysis  is  not  much  more 
complicated  than  for  a  simple  plane  stress  analysis,,  and  it  is 
expected  that  this  new  program  may  find  application  in  future 
investigations  on  other  aspects  of  gravity  dam  behavior* 

The  results  obtained  from  the  crack  closing  analysis  of  the 
Norfork  Dam  do  not  indicate  extremely  severe  stress  values,, 
although  the  maximum  stresses  are  increased  somewhat  over 
corresponding  cases  which  do  not  have  the  full  crack  height. 

In  general,  redistribution  of  stress  which  might  result  from 
cracking  or  creep  effects  in  regions  of  overstress  has  a  favor¬ 
able  influence  on  the  stress  distribution.  Deformation  adjust¬ 
ments  tend  to  distribute  the  stress  more  uniformly  to  the 
entire  system,  and  since  the  structure  is  very  massive,  these 
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average  stresses  are  quite  safe.  Even  the  addition  of  uplift 
pressures  tends  to  have  a  favorable  effect  on  the  system  as 
a  result  of  eliminating  the  region  of  tensile  overstress. 


TABLE  Ilia.  CHECK  OF  BASE  FORCE  EQUILIBRIUM 


Case 

dx 

2* 

3650k 

Error 

dx 

xy 

Fx 

Error 

N-2 

3560k 

-2.5% 

1610 

1530 

5.2 

N~3 

3550 

364C 

-2.5 

1500 

1460 

2.7 

P-2 
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3650 

-0.8 

1650 

1530 

00 

9 

P-4 

3650 

3650 

0 

1530 

1530 

0 

Q-2 

3580 

3650 

-1.9 

1620 

1530 

5.9 

Q-4 

3610 
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-1.1 

1440 

1530 

-5.9 
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APPENDIX 

STRESS  VECTORS  FOR  ALL  CASES 
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